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ABSTRACT 



It is known that the carbon-enhanced, extremely metal-poor (CEMP) stars 
constitute a substantial proportion in the extremely metal-poor (EMP) stars of 
the Galactic Halo, by far larger than CH stars in Population II stars. We in- 
vestigate their origin with taking into account an additional evolutionary path 
to the surface carbon-enrichment, triggered by hydrogen engulfment by the he- 
lium flash convection, in EMP stars of [Fe/H] < —2.5. This process is distinct 
from the third dredge-up operating in more metal-rich stars and also in EMP 
stars. In binary systems of EMP stars, the secondary stars become CEMP stars 
through mass transfer from the primary stars of low and intermediate masses, 
which have developed the surface carbon-enhancement. Our binary scenario can 
predict the variations in the abundances not only for carbon but also for nitrogen 
and s-process elements and reasonably explain the observed properties such as 
the stellar distributions with respect to the carbon abundances, the binary peri- 
ods, and the evolutionary stages. Furthermore, from the observed frequencies of 
CEMP stars with and without s-process element enhancement, we demonstrate 
that the initial mass function of EMP stars need to give the mean mass ~ 10 Mq 
under the reasonable assumptions on the distributions of orbital separations and 
mass ratio of binary components. This also indicates that the currently observed 
EMP stars were exclusively born as the secondary members of binaries, mak- 
ing up ~ 10% remnants of EMP binary systems of mass ~ 10^ Mq in total; in 
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addition to CEMP stars with white dwarf companions, a significant fraction of 
them have experienced supernova explosions of their companions. We discuss the 
imphcations of the present results in relation to the formation of Galactic halo. 

Subject headings: stars: abundances — stars: carbon — Galaxy: halo — stars: 
evolution — stars: luminosity function, mass function 

1. Introduction 

Both HK survey (Beers et al. 1992) and Hamburg/ESO (HES) survey (Christlieb et al. 
2001) have shown that a large fraction (20 ~ 25 %) of extremely metal-poor (EMP) stars 
([Fe/H] < —2.5) exhibit the enhancement of surface carbon abundance. (Rossi et al. 1999; 
Beers 1999; Norris et al. 2001; Christlieb 2003, see also Beers & Christlieb 2005). This forms 
a striking contrast to the fact that more metal-rich CH stars, characterized by enhanced 
CH bands in the spectra, account for only ~ 1% of the other Population II (Pop II) stars 
(Tomkin et al. 1989; Luck & Bond 1991). Figure 1 shows the observed carbon enrichment 
against the metallicity for the stars whose abundances are available from the analyses of high 
dispersion spectroscopy. This figure contains 34 CEMP stars for [Fe/H] < —2.5, and we can 
readily see the prominence of CEMP stars among EMP stars though the sample may not be 
free from bias. In addition, the various peculiar aspects have been revealed; a) CEMP stars 
are divided into two groups, one (CEMP-s) with the enhancement of s-proccss element up 
to 2 ~ 3 dex and the other (CEMP-nos) with normal s-process element abundances (Norris 
et al. 1997a; Aoki et al. 2002a; Ryan et al. 2005, see also Figure 7 below): b) a fairly large 
group of stars display great enrichment of nitrogen relative to carbon (Norris et al. 2002; 
Spite et al. 2005): and c) a bunch of stars show the significant enhancement of r-process 
elements (Sncdcn et al. 1996) with large variations in the degree up to a factor of ~ 100 or 
more (Honda et al. 2004). 

These features may arise from the peculiarities in the structure and evolution of stars 
due to very low metallicity and also from their formation processes, characteristic to the 
early stage of our Galaxy. Theoretically, it has been known that low-mass stars of metal- 
licity [Fe/H] < —2.5 dredge up carbon to the surface by a mechanism different from that 
in the stars of the metal-rich populations (Fujimoto et al. 1990, 2000). This mechanism, 
which is called helium-fiash driven deep mixing (He-FDDM; see §2), is proposed to explain 
the observed features of carbon enhanced EMP (CEMP) stars. Because of this theoretical 
property as well as of the observational characteristics, above mentioned, we call stars of 
[Fe/H] < —2.5 as EMP stars in the following, and distinguish them from "Pop H" stars of 
[Fe/H] > -2.5. 
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There is a well-established scenario for the origin of carbon-rich stars among the Pops I 
and II stars. According to Iben (1975a) and Iben & Renzini (1983), a Pop I or II star enhances 
the carbon abundance in the hydrogen-rich envelope by the third dredge-up (TDU) when 
it evolves to the asymptotic giant branch (AGB). Then the star loses its envelope by the 
stellar wind, eventually turning into a white dwarf. If this star belongs to a binary system, 
the companion star, less massive than this AGB star initially, is exposed to, and accretes, 
the wind matter enriched with carbon and changes its surface composition to exhibit the 
carbon enrichment. As a result, a carbon star remains with an unseen companion star (a 
white dwarf). This is supported by the results of spectroscopic observations that all CH 
stars belong to binary systems with unseen companions (McClure & Woodsworth 1990). 

We propose that all CEMP stars are also made in the binaries similarly to CH stars, 
but with carbon dredged up by He-FDDM as well as by TDU. A statistical analysis of 
carbon-enhanced EMP stars enriched with the s-process elements (CEMP-s) by Lucatello 
et al. (2005a) suggests that most of these stars belong to multiple systems, supporting the 
argument that CEMP-s stars are metal-poor analogs of CH stars. In fact, if plotted as a 
function of the surface carbon enhancement ([C/H]), CEMP dwarfs and CEMP giants have 
distinct distributions, as shown in Figure 2. This is indicative that these stars have suffered 
from the external pollutions and hence changed their surface carbon abundances according 
to the depths of the surface convective zones. CH stars may share a similar feature, but 
because of large metallicity, the giants can hardly display the strong feature of CH hues 
since the dilution of surface carbon by deepening surface convection may reduce the surface 
abundance below C/0 < 1. In fact we see in Fig. 1 that CH giants are predominant in 
the lower metallicity while in the higher metallicity, abundant are Ba stars which are giant 
counterparts of CH stars with C/0 < 1. 

However, some differences are discernible between CEMP stars and CH stars. The 
proportion of CEMP stars to EMP stars is significantly larger than that of CH stars to Pop 
II stars, as stated above. CH stars and CEMP stars have different distributions with respect 
to the orbital separations, as shown in Figure 3, although the statistics may not be sufficient 
for CEMP stars with large separations. In addition, some CEMP stars do not exhibit the 
enrichment of s-process elements, while all the CH stars show the fairly large enrichment. 
The issue of these two subclasses, CEMP-s and CEMP-nos, is raised by Ryan et al. (2005). 
There have been no explanations to where these features come from. 

The first purpose of this paper is to investigate the observational characteristics of 
CEMP stars by compiling data from the literature and to discuss their origins by applying 
the binary scenario, discussed above. The features of CEMP star population have been 
discussed only from the observational viewpoint by Lucatello et al. (2005a) and by Ryan 
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et al. (2005). Relying upon the current understandings of stellar evolution, we present an 
interpretation on the origins of the both subclasses, CEMP-s and CEMP-nos and also of 
the stars enriched with nitrogen but with little carbon. These subclasses are shown to be 
divided by their mass of primary stars. 

Secondly, we examine the binary scenario from the viewpoint of the evolution of binary 
systems with the mass transfer taken into account. We corroborate our interpretation by 
confronting the theoretical predictions with the observed properties of CEMP stars and by 
analyzing the metallicity dependencies through the comparison with CH stars; we deal with 
the relative frequencies, the correlations of orbital period and carbon-enrichment, and the 
distribution of evolutionary stages. 

Finally, we inquire into the initial mass function (IMF) of EMP binaries on the basis of 
our binary scenario. The observed properties of EMP stars are thought to give information 
on the formation history and early evolution of Galaxy. From the viewpoint of chemical 
evolution, the richness of CEMP stars are hnked with the IMF peaking at the intermediate 
stellar mass range (e.g., 4 — 8Mq; Abia et al. 2001). On the other hand, Lucatello et al. 
(2005b) discuss the observed frequency of CEMP-s stars from the binary scenario similar 
to ours, arguing an IMF, shifted shghtly to higher mass as compared with those of metal- 
rich populations, such as expressed in a lognormal form peaking at 0.79 M©. These existent 
studies, however, seem not to properly take into account the recent progress in understanding 
the evolution of currently observed EMP stars, and in particular, the modifications of surface 
characteristics during their long lives. Lucatello et al. (2005b) ignore the difTcrcnccs in the 
evolutionary paths to the carbon enhancement across the metallicity of [Fe/H] ~ —2.5, 
theoretically predicted as stated above. Furthermore, they deal only with CEMP-s stars, 
although it is desirable to put all other constituents of EMP stars into perspective as well. 
In this work, we apply the current best knowledge on the evolution of metal-poor stars and 
consider not only CEMP-s stars but also with CEMP-nos stars and other subclasses. Another 
important factor to the binary scenario is the assumption on the mass ratio distribution 
between the component stars, for which a consensus has not been reached observationally 
nor theoretically (e.g., see Mazeh et al. 2003, and references therein). In this paper, we 
assume one of the plausible assumptions, though different from that assumed by Lucatello 
et al. (2005b, ; see §7), and investigate the consequences. 

The initial mass function derived in this work will give a new insight into the nature of 
EMP stars and their roles in early evolution of our Galaxy and the formation of Galactic 
halo. One proven feature of EMP stars is the scarcity as population in comparison with 
stars of more metal-rich populations. Once Bond (1981) has raised an issue of "G-dwarf" 
problem in the Galactic halo on the basis of the fact that the number of stars decreases 
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sharply with decreasing the metaUicity below [Fe/H] < —1.7 and no stars were found below 
[Fe/H] < —2.6 at that time. Thanks to the recent large-scaled HK and HES surveys, the 
known EMP stars have greatly augmented in number, and yet, it is true that they remain by 
far fewer as compared with other constituent. Pop II stars in the Galactic halo. HES survey 
identifies ~ 200 stars of [Fe/H] < —3.0 in the magnitude range of 10 < -B < 17.5 among 
~ 40% of candidates selected from the fields of 8225 square degrees (Beers & Christlieb 
2005). Accordingly, EMP stars of [Fe/H] < —3.0 may occupy only a tiny fraction of stars, 
much fewer as compared with the stars in the globular clusters, ancillary constituents of 
Galactic halo. We may also address this problem. 

This paper is organized as follows: We start with reviewing the mechanisms to enhance 
the carbon abundance in the stellar envelope of EMP stars in §2. Then we formulate our 
binary scenario in the following two sections; the origin of abundance anomalies observed 
among EMP stars is identified in §3 on the basis of the evolutionary models, and the param- 
eters of binary systems is specified in §4 from the cvohition of binary systems taking into 
account of the stellar wind from the primary and the subsequent accretion by the secondary. 
In §5, the predictions of our binary hypothesis are confronted with the observations from 
particular CEMP stars on the period-carbon enhancement relations. In §6, we investigate 
our binary scenario for GEMP-s stars, whose origin is well-defined theoretically and which 
deserve a separate statistical study because of their relative richness as compared with other 
EMP stars. In §7, we deal with the both subclasses of GEMP-s and GEMP-nos stars to 
derive the IMF of EMP stars and to examine its implications. Finally, §8 concludes the 
paper. 



2. Evolutionciry Chciracteristics of EMP Stars 

A basic distinction in the evolution between the EMP stars and the more metal-rich 
(Population I and II) stars is the engulfment of hydrogen by the helium convection during 
the helium flash in the core or in the shell. Gonsequently, the low- and intermediate-mass 
EMP stars can enrich the surface with the nuclear products of hehum burning through the 
helium- flash driven deep mixing (He-FDDM), triggered by the hydrogen mixing in addition 
to the third dredge-up (TDU) similarly in the metal-rich populations. In this section, we 
review the evolution of EMP stars in relation to the carbon enhancement as well as to the 
s-process nucleosynthesis, and present a revision of the general picture of their evolution, 
originally formulated by Fujimoto et al. (2000) as a functions of initial mass and initial 
metaUicity. 
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2.1. Evolutionary Paths to Surface Carbon-Enhancement 

Fujimoto et al. (1990) first show that the hydrogen mixing occurs during the hehum fiash 
in the core of a metal-free (Pop III) star of mass 1 Mq at the RGB tip, which is consequent 
upon a low entropy of the hydrogen burning shell, realized by low CNO abundances. Hollow- 
ell et al. (1990) compute the progress of subsequent He-FDDM on the RGB (He-FDDM-R) 
in detail and demonstrate that it enhances the surface abundance of carbon and nitrogen. A 
number of following works obtain the consistent results of He-FDDM-R (Schlattl et al. 2001, 
2002; Weiss et al. 2004; Picardi et al. 2004). This mechanism works for the stars of very low 
metallicity ([Fe/H] < —4.5; Fujimoto et al. 1995) and of low masses, M < I.IMq for Z — 
(Suda et al. 2004) and M < 1.2Mq for [Fe/H] ~ -5.3 (Suda 2006 in preparation). 

For the stars that do not undergo He-FDDM-R, Fujimoto et al. (2000) show that the 
similar events happen during helium shell flashes in the early stage of thermal pulsating 
AGB phase. He-FDDM on the AGB (He-FDDM-A) occurs in a star with larger metallicity 
([Fe/H] < —2.5) and of larger masses (M < 3Mq). Cassisi et al. (1996) report the hydrogen 
mixing into helium convection during a helium shell flash for a star of M = 0.8 Mq and 
log Z = —10; despite the sufficiently low metallicity, their model star has skipped He-FDDM- 
R presumably owing to the difference in the input physics affecting the evolution leading to 
the main He flash at the RGB tip. The progress of He-FDDM-A is investigated by a few 
other authors (Iwamoto et al. 2004, a 2 Mq star of [Fe/H] = -2.7; Straniero et al. 2004, 
a 1.5 M0 star oi Z — 5 x 10~^) and all of these results are consistent with Fujimoto et al. 
(2000). 

Figure 4 schematically illustrates the progress of He-FDDM during core flash on the 
RGB tip (upper panel) and He-FDDM during shell flashes on the AGB (lower panel). The 
physical processes are quite similar, and yet, they may differ according to the rate of hydrogen 
mixing when the helium convection flrst makes contact with the hydrogen-containing layer 
(Sweigart 1974). If the mixing rate of hydrogen is sufficiently large, the hydrogen engulfment 
causes the split of convective zone and triggers He-FDDM. For the core helium flash, this is 
the case for the flrst flash because the core flash starts with a sufficient strength, as shown 
in Figure 4 (He-FDDM-R). For the hehum shell flashes on the AGB, on the contrary, the 
hydrogen mixing takes place while the shell flashes grow in strength along with the cooling 
of core, and hence, starts with weak one that cannot bring about the splitting of the helium 
convective zone. After repeating the weak hydrogen mixing as it recurs, flnally the shell flash 
grows strong enough to trigger He-FDDM, as shown in Figure 4 (He-FDDM-A). During the 
He-FDDM(-R and -A), the mixed hydrogen is carried down and burns in the middle of helium 
convective zone and ignites a hydrogen shell flash to split the convection into lower and upper 
ones, driven by the helium and hydrogen burning, respectively. In the decay phase of the 
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hydrogen shell flash, the flash convection retreats in the upper zone and the shell occupied 
by the flash convection expands due to the heat deposited during the hydrogen flash and 
during the helium flash. Consequent upon the shell expansion, the surface convective zone 
deepens in mass and penetrates down into the shells, formerly occupied by the hydrogen 
flash convection. As a result, the nuclear products of the helium flashes, processed by the 
hydrogen shell flash are bought up to the surface. He-FDDM enriches the surface not only 
with carbon but also with nitrogen up to N/C < 1. 

In addition to Hc-FDDM, EMP stars undergo the third dredge-up (TDU) as the stars of 
metal-rich populations except for the very low mctallicity. There is a lower mass limit to the 
stars that undergo the third dredge-up; Lattanzio (1986) shows that the TDU operates in a 
Pop II star of mass 1.5 Mq. A trend is recognized that TDU occurs in less massive stars for 
lower metaUicity (Iben 1983), although the efficiency of dredge-up depends on the treatment 
of convection (see Karakas et al. 2002, and references therein). Furthermore, our low-mass 
models may no longer remain metal-poor as far as CN elements arc concerned because of 
precedent He-FDDM. Thus we may well set the lower mass limit for the TDU to occur at 
~ I.SMq, though this assumption has little to do with the following discussion. This implies 
that stars in the middle mass range between ~ 1.5 Mq and ~ 3.5 Mq experience both, first 
He-FDDM and then TDU (Iwamoto et al. 2004; Straniero et al. 2004). It is also known that 
in massive AGB stars, the temperature at the bottom of envelope convection reaches high 
enough to convert carbon into nitrogen and to reduce the C/0 ratio in the envelope even 
below unity (Boothroyd et al. 1993), while the abundance of sum of CN elements increases 
as carbon is dredged up via TDU. We take the lower mass limit to this hot bottom burning 
(HBB) to be ~ 5M0, although it may depend on the metallicity and on the assumption 
about the efficiency of convective heat transport (see, e.g., Ventura et al. 2001). 



2.2. s-Process Nucleosynthesis in EMP Stars 

The hydrogen mixing may also promote the s-process nucleosynthesis in the helium 
convective zone. The engulfed hydrogen is carried inward and captured by ^^C in the middle 

of convective zone and ^^C, thus prodiiccd, is mixed further inward to release neutrons 
through i3C(«,n)i60 reaction (Iwamoto et al. 2004; Suda et al. 2004). Thus, He-FDDM 
brings about the surface enrichment with s-process elements, as illustrated in lower panel 
of Fig. 4. During the weak hydrogen mixing without the splitting of convective zone, the 
synthesized s-process elements are spread all over the helium convective zone, and will be 
involved again in the helium convective zones during the subsequent helium shell flashes. 
Eventually as the shell flash grows strong enough to trigger He-FDDM-A, s-process elements 
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are dredged up to the surface along with carbon and nitrogen. For massive stars of M > 
1.5Mq, a part of s-process elements, stored in the helium convective zone, are dredged up 
during the following TDU although their abundances suffer dilution by matter newly added 
to the helium core due to quiescent hydrogen shell burning, and then, incorporated into the 
flash convection. On the other hand, He-FDDM-R may bring little s-process elements to the 
surface since the convective zone splits for the first hydrogen engulfment, as shown in upper 
panel, though some of ^^C is carried down into the lower helium convective zone before the 
splitting and gives rise to s-process nucleosynthesis. 

There has been proposed the other site of the s-proccss nucleosynthesis with the same 
reaction ^^C{a,ny^O as the neutron source but attendant on TDU in AGB stars; it relies 
upon the alleged formation of thin layer containing ^^C, i.e., the ^^C pocket, in the top 
of helium zone during the dredge-up phase, presumably by semi-convection in the carbon- 
rich helium zone (Iben & Renzini 1982a,b) and/or by convective overshooting across the 
bottom of surface convection. It is found that the ^^C in the ^^C pocket burns during 
the inter-pulse phase and the s-process nucleosynthesis occurs under the radiative condition 
before the subsequent helium shell flash is ignited (Straniero et al. 1995; Gallino et al. 1998). 
This radiative ^^C burning model is argued to work for the metal-rich populations from 
comparison with solar relative abundances of s-process elements. In particular, a tendency is 
proposed toward larger production of heavier elements for lower metallicity since the number 
of neutrons available per seed nucleus may increase with decreasing metallicity (Busso et al. 
1995, for the review see Busso et al. 1999). Prom this radiative ^^C model, therefore, we 
may expect that EMP stars exhibit very large ratios between the heavy to main s-process 
elements, much larger than the stars of the metal-rich generations. In actuality, however, the 
existent observations indicate that this is not necessarily true; rather reported is an opposite 
trend of decreasing the Pb/Ba ratio for smaller metallicity [Fe/H] ^ —2.5 (Aoki et al. 2000, 
sec also Lucatcllo ct al. 2003). Figure 5 shows the ratios of heavy to main s-proccss elements, 
compiled from the literature, as a function of the metallicity. A break in the variation of the 
Pb/Ba ratio with metallicity is discernible near [Fe/H] ~ —2.5; this is particularly evident if 
we note that two stars with the largest Pb/Ba ratios are among the binaries with the three 
shortest periods, and hence, are likely to have been suffered from an extensive mixing during 
the common envelope evolution (Lucatello et al. 2003). 

The above fact may be interpreted as evidence that the efficiency of ^^C pocket decreases 
and/or varies in EMP stars (Ryan et al. 2001, GaUino et al, 2006, private communication) 
or even as evidence that the radiative ^^C burning will not work for the metallicity below 
[Fe/H] < —2.5 (Suda et al. 2004). The overshooting from the bottom of envelope convection 
may be conceivable in some way or other, but we have not yet a reliable theory for supporting 
the formation mechanism of ^^C pocket, and in particular, for predicting the variations with 
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the metallicity. In addition, for the metallicity of [Fe/H] < —2.5, the s-process nucleosyn- 
thesis occurs in the hehum convective zone triggered by hydrogen mixing within the current 
standard framework of stellar structure, in which mixed ^^C is diluted over the entire helium 
convective zone and will not necessarily give such large numbers of neutrons allotted per seed 
nuclei as predicted from the radiative ^^C burning model (Busso et al. 1999). Accordingly we 
may well assume that there exists a critical metallicity near [Fe/H] ~ —2.5 below which the 
radiative ^^C burning is ineffective during TDU. This assumption will be further examined 
from the comparisons with the observations. 

2.3. A General Picture of the Evolution to Carbon-Enhancement 

The above discussion leads to a general picture of carbon and s-process enhancement in 
EMP stars, as illustrated in Figure 6. It is an up-to-date version of the dependences of the 
He-FDDM and TDU on the initial stellar mass and metaUicity, formulated by Fujimoto et 
al. (2000). 

case I: Low-mass stars of [Fe/H] < —4.5 in the mass range of M < 1.1 M0 ior Z — 
and M < 1.2 Mq for [Fe/H] ~ -5.3 (Suda et al. 2004, Suda 2006 in preparation) 
undergo Hc-FDDM-R at the helium core flash. Once He-FDDM-R occurs, the surface 
composition becomes such as [C/H] ~ and N/C ~ 1. These stars may exhibit little 
enhancement of s-process elements. 

case II and II': For stars located in the region of [Fe/H] < —2.5 and M < 3.5 Mq in 
Fig. 6, but excluding the region of case I, Hc-FDDM-A occurs when they evolve into 
the early phase of thermally pulsing AGB (TP- AGE). The surface abundance of carbon 
reaches as large as [C/H] ~ for the stars of the smallest masses, while it decreases 
for the stars of larger masses because of larger envelope mass and also of smaller mass 
in the flash convective zone. He-FDDM-A ceases to occur after the CN abundances 
in the envelope exceeds [C -|- N/H] ~ —2.5, and the subsequent evolution is divided 
by the occurrence of TDU; case II of M < 1.5 Mq that no further mixing takes place 
and case II' of M > 1.5 M0 that TDU follows and brings C synthesized by hehum 
burning to the surface to reduce the N/C ratio. S-process elements, synthesized via 
the convective ^^C burning, are dredged up by He-FDDM-A for the both cases, and in 
addition, by the subsequent TDU for case II'. 

case III: The stars more massive than ~ 3.5 Mq and with the smallest metaUicity will un- 
dergo neither He-FDDM nor TDU, although the latter point is a matter of controversy 
because of the efficiency of overshooting (Suda et al. 2004). 
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case IV and IV': The stars in these region undergo TDU but not He-FDDM. Cases IV 
and IV' are distinguished by the surface enrichment of s-proccss elements, i.e., case IV' 
of lower metallicity without s-process enhancement and case IV of higher metallicity 
with s-process enhancement. The difference is ascribed to the efficiency of the radiative 
^^C burning. TDU dredges up carbon, but it may be converted into nitrogen via 
the hot bottom burning in the envelope for massive stars of M > 5 Mq; thus, the 
nitrogen abundance may vary with the stellar mass from no enhancement to the large 
enhancement even exceeding the carbon abundance and up to the equilibrium value of 
CN (or CNO) cycles (N/C = 30 ~ 16 in the range of T = 5 • 10^ ~ 10^ K). 

In summary, the surface carbon enhancement in EMP stars is brought about by He- 
FDDM for the initial mass M < 1.5 Mq, by both He-FDDM and TDU for 1.5 Mq < M < 
3.5 Mq, and only by TDU for M > 3.5 Mq. One of the distinguishing features of He-FDDM, 
as compared with TDU, is the surface enrichment of nitrogen at the same time with that 
of carbon; for cases I and II, the surface abundance ratio results N/C~l — 1/5 (Fujimoto 
et al. 2000), while for case II', it decreases as carbon is enriched by TDU. For case IV', the 
nitrogen abundance may range from no enhancement to a large value corresponding to the 
equilibrium in CN (or CNO) cycles. As for s-process elements, our scenario predicts the 
surface enrichment due to the convective ^^C burning for cases II and II', and postulates no 
or little enhancement for case IV'. This conjecture predicts the variations in the abundances 
of nitrogen and of s-process elements in CEMP stars as well. 



3. Binary Scenario I — The Origins of CEMP and Related EMP Stars 

If EMP stars belong to a close binary system, the surface characteristic abundances that 
the evolved primary star has developed can be imprinted onto the secondary star through 
the mass transfer. A CEMP star, thus produced, may keep a record of the primary star on 
their surface characteristics. In this section, we survey the observed characteristics of CEMP 
stars and other related EMP stars, and discuss their origins, in the light of the modifications 
of surface abundances of carbon, nitrogen and of s-process elements, as described in the 
preceding section. 

We compile the characteristics of CEMP and CH stars from the literature as summarized 
in tables 1 and 2, respectively. Our criterion of CEMP stars is that [Fe/H] < —2.5 and 
[C/Fe] > 0.5, but we include those with shghtly larger metaUicity which show the nitrogen 
enhancement ([N/Fe] > 0.5 ) with taking into account possible errors (~ 0.25 dex) in the 
abundance determination. We have 34 CEMP stars in total (see Fig. 1), fist their available 
data of the orbital parameter, the gravity, the metallicity, the abundances of C, N, O and 
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s-process elements (Sr, Ba, Pb) relative to iron; 7 stars have the orbital periods derived with 
the maximum P = 11.3 yr and 7 more stars are suspected of binarity from variations of 
radial velocity. As for CH stars, we compile 85 stars by the criterion that [Fe/H] > —2.5; in 
addition to those classified as CH stars in literature, we include three blue stragglers showing 
strong CH lines, and also, Ba stars which are thought as giant counterparts of subgiant CH 
stars with C/0 < 1 (e.g.. Luck & Bond 1991). 38 stars of them have the orbital periods 
determined (two stars with the lower bounds). Note that CH stars are divided into two 
groups of CH (giant) stars and subgiantCH stars, the latter of which includes both dwarfs 
and subgiants (Bond 1974). 



3.1. Sub-classification of CEMP Stars 

Figure 7 shows the observed relationship of the enrichment between Ba and nitrogen 
among CEMP stars. As already noticed by Aoki et al. (2002a, see also Ryan et al. 2005), 
they are rather clearly divided into two groups of [Ba/Fe] > 1 and [Ba/Fe] < 0, where a few 
exceptions include CS22892-052 enriched with r-process elements (Sneden et al. 1996). 

All CEMP stars with Ba enrichment ([Ba/Fe] > 1), which are defined as CEMP-s, are 
enriched with nitrogen ([N/Fe] > 1). The relative abundance to carbon lies in the range of 
0.01 N/C < 1, and is consistent with the theoretical predictions of cases II and II'; the 
largest ratios may result from He-FDDM in case II and the smaller ratios are exphcable in 
terms of the increase of carbon via TDU in case II'. 

On the other hand, CEMP stars without Ba enrichment ([Ba/Fe] < 0), which are 
defined as CEMP-nos, exhibit still larger variations in the nitrogen abundances, ranging 
from no enhancement [N/Fe] ~ (CS 22877-001) to much larger enhancement with the 
N/C ratios ~ 10 (CS22949-037). The range of enrichment well fits to case IV' evolution in 
our scenario with and without hot bottom burning. It is difficult to search for the site of 
nucleosynthesis in the stars other than hot bottom burning that can realize such large N/C 
ratios as near to the equilibrium values in the CN cycles. As a class, CEMP-nos stars have 
the carbon enhancement smaller than CEMP-s stars, as pointed out by Aoki et al. (2002a, 
also see Fig. 7). This is attributable to the dilution of carbon, dredged up, in the envelope 
of the primary star, which is more massive in case IV' than in cases II and II', and also, to 
the transformation of carbon into nitrogen by hot bottom burning. 

Some of stars in the CEMP-nos subclass have large oxygen abundance, comparable to, 
or even larger than, the carbon abundance. In the binary scenario, the enhancement of 
oxygen exceeding about a tenth of carbon abundance has to be attributed to the pristine 
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abundance (cf. Suda et al. 2004), which may well be explained in terms of the general 
tendency of [0/Fe] increasing with decreasing metallicity and its extension to EMP stars 
(Israelian et al. 1998; Boesgaard et al. 1999, but see also Garcia Perez et al. 2005) (see 
appendix A.l). 

Among EMP stars, a class of stars display large nitrogen enrichment with [C+N/Fe] > 
but display only weak or little carbon enhancement with the ratio C/0 < 1, which are 
classified as "mixed" stars by Spite et al. (2005). We propose that these stars are made 
by the same mechanism as CEMP-nos stars but carbon was converted into nitrogen by hot 
bottom burning in more massive primary stars than that of CEMP-nos stars (appendix A. 2). 
In our compilation, those with r-process clement enhancement such as CS22892-052 may also 
be included in the subclass of CEMP-nos since the Ba enrichment can be attributed to the 
r-process nucleosynthesis (Sneden et al. 1996). We will inquire their origin separately in this 
paper (§7.3 and appendix A. 3). 

In our binary scenario, the origins of subclasses of EMP stars are divided according to 
the mass of their primary stars. CEMP stars enriched with s-process elements (CEMP-s) 
are formed in the binary system with a primary star of mass in cases II and II' while CEMP 
stars without the s-process element enrichment (CEMP-nos) in the binary system with a 
primary of mass in case IV'. The nitrogen- enhanced "mixed" stars may be assigned to the 
binary systems, same as CEMP-nos stars, but with more massive primaries. 

4. Binary Scenario II — Binary Parameters of CEMP Stars 

In this section, we discuss the binary parameters that the secondary star can evolve to 
CEMP stars through the mass transfer. The first condition is that the primary star is allowed 
to develop the surface carbon enhancement. The second one is that the binary system can 
enrich the secondary star sufficiently with the envelope mass ejected from the primary stars. 
These conditions set lower and upper bounds to the binary separations, respectively, for a 
given set of the masses of primary and secondary stars, as discussed below. We can derive 
the relation of the carbon abundance of secondary star to the mass of primary stars and the 
orbital separations (periods), which may be compared with the observations. 

4.1. Roche Lobe Overflow 

We start with the discussion of the first condition. In close binary systems, the primary 
star expands as it evolves, and sometimes, fills its Roche lobe to lose the envelope through 
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the mass transfer to the secondary and/or the mass out- flow from the systems (e.g., Iben 
& Tutukov 1985). If this Roche lobe overflow (RLO) happens before the star becomes a 
carbon-enhanced star, then this system cannot yield a CEMP or CH star. In order to yield 
a CEMP-s star, therefore, the initial size of the Roche lobe has to be greater than the stellar 
size when the primary star for M < 3.5 Mq starts the He-FDDM. The same is true for a 
CEMP-nos star and for a CH star but with the stellar size when TDU starts. Since the 
volume of Roche lobe is given as a function of the mass ratio M2/M1 and the separation A 
(Paczynski 1971), the lower bounds, ^He-FDDM and Atduj to the initial binary separation 
necessary to accommodate the primary stars when He-FDDM and TDU start, may be derived 
from the radii, -RRe-FODM and -Rtdu, at the stage of He-FDDM and TDU as; 

. / A \ -RHe-FDDM(or i?TDu) x 

>lHe-FDDM(orATDu) - 0.38 + 0.2 log(Mi/M,) " 

Figure 8 shows these lower bounds on the diagram of the initial stellar mass of primary star 
and the binary separation for the stellar models taken from existent computations (Iben 
1975a; Lattanzio 1986; Marigo 2002; Suda et al. 2004, , Suda 2006 in preparation). The 
initial orbital radius needs to be greater than 0.2 ~ 1 AU for a CEMP-s star while greater 
than ~ 1 AU is necessary for a CH star and for a CEMP-nos star. 

After the primary has developed the surface carbon-enrichment, the RLO may occur 
on AGB before the primary evolves to a white dwarf. In this case, the binary will undergo 
common envelope evolution and shrink in size to be a short-period system or even to coalesce 
in some cases. These systems may be distinguished from CH and CEMP stars produced 
through the accretion process in the wind from the primary star in the binary of wider 
initial separation as proposed by Boffin & Jorissen (1988) and discussed in the following. 



4.2. Wind Accretion from AGB Stars 

The second condition, i.e., an upper bound to the binary separation, results from the 
wind accretion. AGB stars in binary lose their envelope through stellar wind (Reimers 
1975; Gail & Sedlmayr 1987) if the separation are sufficiently large to avoid RLOs. Then, 
the secondary stars may accrete carbon-rich material in this wind to be carbon-enhanced 
stars. The process for the secondary stars to accumulate the gas from the wind may well 
be approximated to the Bondi accretion (Bondi & Hoyle 1944), and the cross section, aacc, 
may be evaluated from the following formula. 



(2) 
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where G is the gravitational constant, M2 the mass of the secondary star, and Vre\ the velocity 
of the wind relative to the secondary. Thus the mass accretion rate dM2/dt of the secondary 
star is related to the mass loss rate dMi/dt of the primary as 



dM2it) 



(3) 



dt A{t)^Vrei{t)^ v^nd ' ' dt 

where the separation A(t) between the two stars has been introduced as a function of time, 
t, and t'wind is the wind velocity with respect to the primary star. The relative speed of the 
secondary star to the stellar wind is given by 



wind 



+ Vl 



V2 



(4) 



with the orbital velocities, vi and V2 of primary and secondary stars. The mass accreted by 
the secondary star depends on the separation A{t) and the relative velocity Vrei{t), which 
may vary during the mass transfer along with the masses Mi{t) and M2{t). 

Suppose that the primary and secondary stars rotate on each circular orbit around 
the center of mass while the primary loses mass at a constant rate dMi/dt, then the basic 
equations governing the evolution of a binary system are 

M,{t)v,{tf/n{t) = M2{t)vl{t)/r2{t) 

= GMi{t)M2{t)/A\t), (5) 

Mi(t)ri(t)-M2(t)r2(t) =0, (6) 

n{t)+r2{t) = A{t). (7) 

Here denotes the radius of the orbit of star i and Vi the velocity. The subscript i = 1 refers 
to the values of the primary star and i = 2 to those of the secondary star. The timescale for 
the change of A{t) is assumed to be much longer than the orbital period. 

If Jeans' theorem is applied, the relation between the separation and the mass of the 
binary system is given by 

A(t)[M,{t)+M2(t)]^X, (8) 

where A is a constant, given by the initial conditions at time tinitiai- This implies that the 
wind is assumed to carry the angular momentum away from the system without transferring 
it to the other material. By expressing A{t) and Vrei(^) in terms of the stellar masses and 
the wind velocity, equation (3) is reduced to 



dM2it) 
dMi{t) 



M2{t) 



Mi{t) + M2{t) 



X 



1 + 



M^{t) + M2{t) 
-3/2 



{M,{t) + M2{t)y 



(9) 
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where is a constant, defined as 

= \[><JG = V2Mi(Mi + M2)(yl/i?i) Kind/^^esc] ■ (10) 

In the most rightmost member, we normahze the wind velocity with respect to the escape 
velocity, Vesc, from the surface of the primary star (= ^j2GMi/ Ri, where Ri is the radius of 
primary star). 

If A » then, Ijl/{Mi + M2) » 1 (since t'wind ~ ^esc) so that the variation of M2 
is much smaller than that of Mi. We may neglect the variation of M2 in the right-hand 
side and integrate eq. (9) until the primary becomes a white dwarf at time ignai to find the 
accreted mass. Mace, on the secondary star; 



Mace 



M| / Ml + M2 , Y 1 



-log - (Ml + M2 



+ V/^2 + (Mi+M2)2))) (11) 
/// JMl=Ml(ti„itial). 

The relation between the initial stellar mass, Mi(iinitiai), and the mass, Mi(ifinai), of C-|- O 
core for metal-poor stars may well be approximated to 

Mi(tfinai) = max[0.54 + 0.073Mi(iinitiai), niin[0.29 

+0.178Mi(tMtiai), 0.65 + 0.062Mi(ti„itiai)]], (12) 
that for the metallicity of Zq — 0.001, given by Han et al. (1995). 



4.3. Dilution due to Surface Convection 

Suppose that the matter of constant carbon abundance, -'^cp, is accreted onto the sec- 
ondary of the initial carbon abundance Xc,2,o- If the mass in the surface convection, Afscz.D, 
remains constant through mass accretion, then, the variation of surface carbon abundance, 
Xc,2, of the secondary star may be described by; 

dXc,2 ^ Xc,p - Xc,2 dM2 

dt Msez,D dt ' 

The integration yields Xc,2 as a function of the accreted mass, Maee, as 

Xc,2 = Xc,2,0 + (^C,p - Xc,2,0) [1 - exp(-Maee/M,e.,D)] • (14) 

The mass transfer is likely to occur while the secondary star resides on the main sequence 
with a shallow surface convective zone. Then, the surface abundance suffers an additional 
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dilution owing to the deepening of surface convection as the secondary later evolves to the red 
giant branch. For the metallicity of [Fe/H] < —3, a star of mass M2 — 0.8 Mq has the surface 
convection of mass Mscz,d — 0.003 Mq on the main sequence. The surface convection, though 
once decreasing as the hydrogen abundance decreases in the center, turns to deepen rapidly 
as the shell burning starts and attains at the deepest reach of Mscz,g — 0.20 ~ 0.35 M0 in 
mass near the base of red giant branch when the luminosity L ~ several x 10 Lq (Fujimoto 
et al. 1995; Suda et al. 2004). At the base of RGB, therefore, the surface carbon abundance 
reduces to be 



(here we assume Msc7,,g > Mice + Afscz.o)- After that, it remains constant along the giant 
branch. As a corollary, if carbon is brought by external pollution, the [C/H] value decreases 
up to the maximum of 2 dex when an EMP star evolves to a giant. 

The surface carbon abundance of secondary stars resulting from the wind accretion is 
given as a function of the initial masses of component stars and the initial binary separation 
from eq. (11) and from eq. (14) or eq. (15) for dwarfs and giants, respectively. These 
relations may be transferred to the upper bound, Am(Mi, M2, [C/H], Mgc^), to the initial 
binary separation that allows the binary system with the primary star of mass Mi to enrich 
the secondary star of mass M2 that has the surface convection of mass M^cz with carbon 
above [C/H] through the wind accretion. In Fig. 8, we plot the upper bounds of initial 
binary separation that give the carbon enhancement of [C/H] = —1, —2, and —3 for the 
secondary of giants with the surface convection Mgcz = 0.35 Mq (M2 = 0.8 Mq). The carbon 
enhancement necessary for CEMP and CH stars may be taken to be [C/H] > —3 and — 1, 
respectively. For such a large initial separation, since the orbital velocity is much smaller 
than the wind velocity, the cross section of Bondi-Hoyle accretion reduce to a constant, 
(Tacc = 7r(G'M2/i'wind)^) ^-nd hcncc, we may approximate 



The upper bound of separation increases in inversely proportion nearly to the second power 
of the carbon enrichment of secondary star, and hence, CEMP stars may have much wider 
separations than CH stars. The bound also depends on the carbon abundance, [C/H]p, in 
the wind. For TDU, the latter is currently subject to uncertainties for lack of the reliable 
theory on the mass loss, and yet, we may well assume that [C/H]p ~ or so. 

Now we have formulated our binary scenario as a function of the mass ratios and the 
binary separations, as illustrated in Fig. 8. In the following three sections, we examine the 
validity of our interpretation and binary scenario firstly as to the relationship between the 



Xc,2 = ^C,2,o(l - Ma^c/Mscz.c) + Xc,^{M^/ M^^^q) 



(15) 



Am(Mi, [C/H], Msez) = [{Ml - Mi(tfi„,i)} 




(16) 
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orbital period and carbon enhancement, secondary as to the statistical properties of CEMP-s 
stars, and finally, as to the relative frequencies of CEMP-s and CEMP-nos stars and discuss 
the initial mass function of EMP stars. 



We can deduce the relations between the carbon abundance of the secondary star and 
the orbital period from our modeling in the preceding section. The accreted mass, Mg^c, 
onto the secondary star can be estimated from eq. (11) as a function of the initial masses of 
component stars, Mi(tinitiai) and M2(tinitiai) and the initial separation, A{tinitia.i) , and then, 
the surface chemical compositions of the secondary stars may be evaluated with the mixing 
in the surface convective zone of the secondary star taken into account from eq. (14) or 
eq. (15). On the other hand, as a kinematic feature, the orbital period after the primary 
star becomes a white dwarf is calculated from the final separation A(tfinai) in eq. (8) and the 
final mass Mi(ifinai) of the primary star in eq. (12). These are observable and can be used 
to test the present scenario, as in the similar way to Boffin & Zacs (1994) who discuss the 
relationship between the overabundances of s-process elements and the orbital periods for 
barium stars. 

In deriving the relationship between the accreted mass and the orbital period, we have 
applied Jeans' theorem. In actuality, however, the wind material is thought to carry an 
additional orbital angular momentum unless the wind velocity Vwind is much larger than the 
orbital velocity forb (= l^^i — ^21). In fact, hydrodynamical calculations of wind accretion 
(e.g., Nagae et al. 2004; Jahanara et al. 2005) have shown that the orbit shrinks if Wwind ^ 
(1.5 ~ 1.7)forb ioT the mass ratio of q{= M2/M1) = 1 ~ 1/3 (see also Hachisu et al. 1999), 
which may correspond in our scenario to; 



Jeans' theorem holds good for the wider binaries of Vwind ^ (1-8 ~ 1.9)i;orb; while the flow 
structure resembles the Roche lobe over-flow (RLO) for the wind velocity of Vwind < 0.4 Vorh- 

Figure 9 plots the relations between the surface carbon abundance [C/H] of secondary 
stars and the orbital period, calculated under the assumption of Jeans' theorem with the 
primary stars of mass 1.5 and 3 Mq. The bottom panel is for the binary systems of [Fe/H] > 
—2.5 whose primary stars have undergone the carbon enhancement by TDU, and may be 
compared with the observations for CH stars. Here we assume the carbon abundance in the 
wind matter to be [C/H] p = 0.5 and the wind velocity to be I'wind — 20 kms~^. A subgiant 



5. 



Period-Carbon Enhancement Relation 



A < 5 ~ 7 (1 + g) (Ml/Mo) (T;^ind/20 km s-^)-^ (AU). 
P < ll~16(l + g)(Mi/Mo)(^;^i„d/20kms-i)-3(yr). 
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has thinner surface convective zone (by a factor of 30 in mass) than a red-giant star, which 
yields larger carbon enhancement with the same binary parameters. For the systems of 
short period -P ^ a few x 10 yr, the subgiants tend to have the same surface abundances 
as the wind. On the other hand, the giants display the carbon abundance smaller than 
the wind even for the smallest separations, and decreases for larger separations nearly as 
[C/Fe] ~ —4/3 log P. These tendencies agree with the observations that the subgiants tend 
to have larger carbon abundance than the giants; the exceptions are the blue stragglers of 
smallest metaUicity of [Fe/H] < —2, which may be attributed to the difference in the carbon 
abundances reached in the primary stars, as discussed below. For most of subgiant CH stars, 
the observed carbon abundances are larger than the solar value, which indicates that the 
carbon abundance in the wind is above the solar value (up to [C/H]p ~ 0.5). The observed 
scatter may be taken as reflecting the different efficiencies of TDU to bring carbon into the 
envelopes of primary stars and/or attributed to the gravitational settling of accreted carbon 
in the envelope of secondary stars. Weiss et al. (2000) show that the heavy pollutants settle 
down to decrease their surface abundance by a factor up to ~ 10 in 10 Gyr, which may 
partly explain the spread of surface carbon abundances, observed from CH sub-dwarfs and 
also CEMP stars in Fig. 2. The gravitational setthng has httle thing to do with the carbon 
enrichment of giants, however, since it is overtaken by the surface convection deepening along 
the red giant branch. But note that for [Fe/H] > —6, the surface convection reaches deep 
enough to the shell where the hydrogen shell-burning passes, and hence, the accreted matter 
may affect the evolution of the secondary to the helium core flash (Fujimoto et al. 1995). 

As for the separation (or the period), the observations hes well below the predictions 
based on Jeans' theorem. Since the periods of the observed systems most fall in the range 
of eq. (18) or shorter, they are thought to have suffered the orbital shrinkage during the 
mass transfer episode. In order to see the effects of additional angular momentum loss, we 
flgurc out and plot the relations under the assumption of constant separation in the flgure 
(denoted by broken lines), too; they end with smaller flnal separation and larger accreted 
mass as compared to the relations based on Jeans' theorem. The observed points mostly 
lie near the shortest ends of the separation necessary to contain the primary stars within 
the Roche robe at the beginning of TDU. Furthermore, some of them including the two 
blue stragglers have the observed periods even shorter than the latter, which suggests that 
these systems have experienced a large shrinkage of separations during the common envelope 
phase. 

The top panel shows the results for CEMP stars, where the carbon abundance of the 
wind is taken to be [C/H]p = 0. Since the radius of EMP primary at He-FDDM is smaller 
than the radus of Pop 11 primary at TDU (see Fig. 8), the binary separation can be smaller 
for CEMP stars than CH stars. CEMP stars of short period P < 10 yr must also have 
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suffered the shrinkage of separations during the mass transfer episode. If the shrinkage is 
taken into account, the relation calculated from the wind accretion scenario gives reasonable 
description to the observed trend of CEMP stars similarly to CH stars, discussed above. It 
is worth noting that CEMP-nos star, CS22957-027 ([C/H] = -7.1 and P = 3,125 days) 
is also among the stars showing reasonable agreement with the theoretical prediction. A 
dwarf HE0024-2523 of the shortest period, and possibly subgiant star CS29497-030 (names 
attached in the figure) , have undergone RLO to accrete carbon-rich matter from the primary 
stars. The eccentricities of these stars are nearly zero, which is also expected from RLO, 
though the detailed modeling of their evolution awaits future studies with the interaction 
between the stellar wind and the orbital motion properly taken into account. The other dwarf 
G77-61 of the second shortest period may deserve comment because of rather weak carbon 
enrichment of [C/H] ~ —1.4 (Plez & Cohen 2005). It has the mass ~ 0.3 Mq (Dearborn et 
al. 1986) and is likely to be wholly convective. In addition, the cross section of the star for 
the accretion is a factor of ~ 7 smaller than that of a red- giant CEMP star [see Eq. (2)], 
and hence, the accreted mass of carbon-rich material might be smaller by the same factor. 
Even if we take into account these two factors that may lead to a relatively low value of 
[C/H], our model is difficult to reproduce the values of G77-61 with the carbon abundance 
[C/H]p = 0. Accordingly, this star requires both the shrinkage of orbital separation and 
the carbon abundance in the wind from the primary smaller than the solar abundance 



One of the factors that may affect these correlations is the difference in the CNO abun- 
dance in the wind between CEMP and CH stars, as discussed above. The difference between 
CEMP dwarfs and CH dwarfs as large as ~ 0.5 dex is discernible also from the distribution 
of carbon abundances in Figs. 2. This may reflect the difference in the efficiency of TDU 
due to the metallicity; possible reason may be the increase in the hydrogen shell burning 
rate and the resultant expansion of envelope with increasing surface enrichment of carbon. 
Another factor may be the eccentricity of the orbit. The effect of the orbital eccentricity 
on the mass accretion rate can be incorporated as was done by Han et al. (1995), and a 
larger eccentricity tends to increase the mass accretion rate and thus increase the carbon 
abundance of the secondary. 



([C/H]p 



1)- 



6. 



Statistical Features of CEMP-s Stars 



We have proved that our binary scenario gives a reasonable description of the observed 
characteristics for the particular CEMP stars with known binary periods in comparisons 
with CH stars. Now we turn to the statistics of CEMP stars to demonstrate the vahdity of 
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our binary scenario and to investigate insights into their properties, gained from our binary 
scenario, in the following two sections. We first deal with CEMP-s stars of case II and II', 
which are well delineated as a subclass, both observationally and theoretically by the large 
enrichments not only of carbon and nitrogen but also of the s-process elements. We then 
demonstrate that constraints can be drawn upon the initial mass function and the nature of 
EMP binaries with taking both CEMP-s and CEMP-nos stars into account. 



On the basis of our scenario, we may estimate the fractions that CH and CEMP-s stars 
occupy among the Pop II and EMP stars, respectively. The binary parameters that the 
secondary stars can evolve to CEMP-s and CH stars arc depicted on the diagram of initial 
mass of primary stars and binary separation in Fig. 8. For CEMP-s stars, we may take 
the lower bound to the initial mass of primary stars to be Mi(tinitiai) = 0.8 and assume 
that the primary stars more massive than this mass have already ejected their carbon-rich 
envelope and evolved to white dwarfs. The upper bound to the primary initial mass for 
CEMP-s stars is set at MHe-FDDM — 3.5 Mq by the condition for He-FDDM to occur. For 
CH stars, on the other hand, the lower and upper bounds are set by the lower mass limit 
to TDU at Mtdu — 1-5 Mq, and by the hot bottom burning (HBB), which converts carbon 
into nitrogen in the surface convection, at Mhbb — 5Mq, respectively. As for the initial 
binary separation, the first condition is that the primary star is allowed to develop carbon 
enhancement and sets the lower bounds, ylne-FDDM and Atdu, for CEMP-s and CH stars, 
respectively. The second condition is to enrich enough the carbon abundance of the secondary 
star through the wind accretion, which imposes the upper bounds, Am{Mi, [C/H], M^cz)- The 
carbon enhancement necessary for CEMP and CH stars may be taken to be [C/H] > —3 and 
— 1, respectively. Consequently, CEMP-s stars occupy much larger parameter space than 
CH stars, as seen from Fig. 8. 

In order to evaluate the relative frequencies of CEMP-s and CH stars numerically, we 
have to specify the initial mass function ^(Mi) for the primary stars and the distributions, 
n{q) and f{P), of binary systems with respect to the mass ratio q (= M2/M1) and the orbital 
period P, respectively. The fractions of CEMP-s and CH stars to EMP and Pop II stars, 
0CEMP-S and 4>cR, are then proportional to the following integrals: 



6.1. Frequency in Comparison with CH Stars 
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/-Mhbb ^(M2/m) ^ 

./mtdu ^ 

fAM(m,[C/H],Mscz) 
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As for the initial mass function (IMF) for the primary stars of mass > 0.8 M©, we may here 
assume Salpeter's power law, 

e(m) oc (21) 

We may also assume a lognormal form as in the next section (see blow eq. (36). For the 
binary periods, we may adopt the distribution which has been derived from the observations 
of local stars, 

(Duquennoy & Mayor 1991), where P is the period in units of days. 

On the contrary, as for the mass ratio distribution, no consensus has been achieved. For 
the near-by G-dwarf binary samples, it is shown that the distribution function can well be 
fitted by a Gaussian form, centered at g = 0.23 with a rather large dispersion of ~ 0.42, or 
could be even flat for q < 0.23. It is also pointed out, however, that the same mass function 
as for the single stars in Millar and Scalo law is admissible for the secondary stars because 
of large uncertainty in the mass ratio below q < 0.1 (Duquennoy & Mayor 1991). Even with 
the recent studies it still remains an open question. On one hand it is reported that the 
distribution have a small peak at q ~ 0.8 and rises as the mass ratio decreases to q ~ 0.2 
with a drop at g < 0.2 for nearby spectroscopic binaries (Goldberg et al. 2003) or that it 
has a broad shallow peak between q ~ 0.2 and ~ 0.7 and a sharp peak at g > 0.8 for nearby 
F7-K stars (Halbwachs et al. 2003). On the other hand Mazeh et al. (2003) argue that the 
mass ratio distribution approximate constant over the range g = 1.0 — 0.3 and a constant 
distribution cannot be ruled out at lower g for main-sequence binaries measured by infrared 
spectroscopy. In particular, Goldberg et al. (2003) find that the halo binaries have a flatter 
distribution than the disc binaries. For such low mass binaries, the mass-ratio distribution 
is subject to large uncertainties in the range of mass ratio below g ~ 0.1 since the main 
sequence does not exist for M < 0.08 Mq. For massive star binaries with 0-type and B-type 
primaries, in contrast, the mass function of the secondary stars is found to extend into the 
very small mass ratio g < 0.1 (Abt 1983). In the present study, therefore, we may well 
assume a uniform distribution of mass ratio of binary systems for simplicity; 

n(g) = const for < g < 1.- (23) 

The resultant mass functions of CEMP and CH stars depend only weakly on the secondary 
mass M2 through the upper and lower boundary values of integral with respect to the 
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separation. This makes a sharp contrast with the case when the same single-star mass 
function as the primary stars is assumed for the secondary stars, as done by Lucatello et 
al. (2005b). We may examine the vahdity of our assumption from the comparisons with the 
observed properties of CEMP and CH stars, as discussed later. 

The observed flux-limited samples of CEMP-s stars and CH stars are dominated by 

those which have masses at nearly upper end and are most luminous. Putting M2 — 0.8 Mq 
and evaluating the double integrals with respect to the period and primary mass in eqs. (19) 
and (20) numerically with [C/H] = —3 and —1 respectively for CEMP-s and CH stars, we 
find that: 

0CEMP-s/0CH ^ 14. (24) 

This number represents the ratio between the frequencies of CEMP-s and CH stars among 
the stars of corresponding metallicity that are born in binaries. Therefore, it gives reasonable 
account of the observed much larger frequency of CEMP stars among EMP stars as compared 
with the frequency of CH stars among stars of Pop H stars, unless the fraction of EMP stars 
born in binaries is significantly smaller than those of the more metal-rich populations. In 
evaluating the ratio, we have apphed the Salpeter's IMF to EMP stars as well as to CH stars, 
although there is no guarantee that both obey the same initial mass function. Nevertheless, 
this conclusion will be little affected since the large ratio arises essentially from the difference 
in the parameter space that EMP and Pop 11 stars can develop the carbon-enhancement, as 
seen from Fig. 8. We will discuss the influence of different IMFs in the next section. 



6.2. Period Distributions in Comparison with CH Steirs 

Since CEMP-s stars and CH stars stem from the binary systems of different parameters 
in our binary scenario, they should differ in the distribution of the orbital separation (or the 
period). The relative frequencies of the initial orbital separation (or the initial orbital period) 
may also be predicted from the integration of eqs. (19) and (20) with respect to the primary 
mass for CEMP stars and CH stars, respectively. The derived relative frequencies are plotted 
in Fig. 3 and compared with the observations. Here, we take into account both giants and 
dwarfs and a shoulder on the right-hand side of each curve comes from the contribution of 
dwarfs. 

For CH stars, our binary scenario predicts the initial separations (or periods) restricted 
to rather narrow range between 3 ~ 30 AU (several ~ 100 years). The observed periods also 
fall in a narrow range of separations (or orbital periods) in agreement with the prediction, 
but are shifted toward smaller separations (or shorter periods) by several factors on average. 
This shrinkage of binary systems is consistent with the argument in the preceding section. 
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based on the eqs. (18), which anticipates that the wind may carry away larger amount of 
angular momentum than specific to the orbital motion of primary star in interaction with 
the secondary motion. 

In contrast, CEMP-s stars are predicted to stem from the binaries of much larger range 
of initial separations (or initial periods); for those with [C/H] > —3, the initial parameters 
range from 0.3 AU (less than a year) to more than 1000 AU (~ 10000 years) with the central 
values near ~ 30 AU (~ 100 years). All the observed orbital periods of CEMP-s stars are less 
than ~ 10 yrs because of limited span of monitor time to date, and fall within the shorter 
half of initial period range, predicted from our scenario. For such shorter periods, the binary 
systems must have suffered from shrinkage during the mass loss event from the primaries. 
The observations are consistent with this expectation, where the distribution of observed 
separations (or periods) is extended toward shorter separations (or periods) than predicted. 
In addition, we expect from Figs. 3 and 8 that more than a half of CEMP-s giants, and still 
larger portion of CEMP-s dwarfs, have wider separations (or longer periods) than observed 
to date. In this parameter range, Jeans' theorem tends to hold and the binary separations 
(or periods) further increase. So, in the extreme end, it is difficult to detect the orbital 
motion of binary spectroscopically as a periodic variation of absorption lines. 

In summary, the binary scenario is compatible not only with the observed period distri- 
bution of CEMP-s stars but also with the lack of detected variations in the radial motion for 
a half of them to date. These results are indifferent to the assumptions on the initial mass 
function and on the distribution of mass ratio. They will not be basically influenced by the 
assumption on the period distribution, f{P), either, unless it varies greatly in the range of 
separation, defineated as the parameter space for CEMP-s stars in Fig. 8. 

6.3. Frequencies of Dweirfs and Giants 

CEMP stars make a sharp contrast with CH stars in their relative distributions with 
respect to the evolutionary stages, too. Among CEMP-s stars, giants ovcrwcigh dwarfs and 
subgiants in number (14 giants among 21 stars in our sample), while subgiant CH stars 
outnumber CH (giant) stars (18 CH giant stars and 25 subgiant CH stars in our sample). 
From the binary scenario, this may be explicable in terms of the effect of the pristine metal 
abundances in the dilution of carbon-enrichment when a star evolves to be a giant. For such 
large pristine metallicity as CH stars, the mixing of accreted carbon-rich matter into the 
extended surface convection of giants is ready to convert the carbon-rich dwarf stars to the 
giant stars exhibiting no enhancement of carbon, i.e., C/0 < 1; we see in Fig. 1 that CH 
giants are predominant at lower metallicity; if Ba stars of C/0 < 1 are included, giants may 
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exceed dwarfs in number (see Fig. 2) although we should take into account the differences 
in the metallicity dependence (e.g., sec Smith & Lambert 1986). For CEMP stars, on the 
contrary, the carbon enrichment may be discernible even if drcaded-up carbon is diluted in 
the surface convection of giants of mass larger by a factor of ~ 100 than the dwarfs. 

In the binary scenario, the mass distribution of CEMP-s stars is expressed from eq. (19) 

as; 



;^He-FDDM n(M2/m) 

^CEMP-s(M2) = / dmCim) ^ " ' 



0.8 Mq 



m 



/■Am (m, -3, Mscz) 

X / f{P)-rda. (25) 

Wc may estimate the ratio between giants and dwarfs from our scenario by using the results 
of evolutionary study of a star of mass of 0.8 AIq for the metallicity [Fe/H] = —3 (Suda 
2006 in preparation). It takes U = 15.9 Gyr for this star to evolve from the zero-age main 
sequence to a white dwarf, including 1.3 Gyr (= Atsc = tg — ^to) spent as subgiant after 
the turn-off (at t — ^to) to the stage when the surface convective zone becomes deep (at 
t — tg), and then, 0.45 Gyr (= Ato — t{ — tg) before it becomes a white dwarf at t — tf. The 
mass ranges, Amso and Amc, of subgiant and giant stars that were born in early days (~ tf 
years ago) can be estimated from the mass-luminosity relation for dwarf stars and the time 
spans during which the star resides as a subgiant and on the RGB and horizontal branch, 
respectively. The mass-luminosity relation for dwarf stars is assumed to follow a power law 
-^dwarf('^) oc (ct = 3.5). From the relation 

Ata/tf = (a - 1) (Am/O.SMo) , (26) 

we obtain Amc = 0.01 Mq and Amsc = 0.03 M©. 

Since the observed EMP stars constitute flux-limited samples, we should also take into 
account the differences in the luminosity among them and their spatial distributions. We 
denote, by Ns{L), the surface number density (per unit steradian) of stars with the luminosity 
L that are accessible with a given limiting flux /umit- As for the special distribution, we may 
well assume that the EMP stars obey the density distribution, p(r) , of stars in the Galactic 
halo, approximated by 

p(r) oc r-^ (27) 

as a function of radial distance r from the Galactic center (e.g. see Majewski 1993, and 
references therein). Since the surveys have been performed exclusively in the flelds of high 
galactic latitudes, we may well consider the direction perpendicular to the Galactic disk. 
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which leads to; 
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(L) ^ / pir)z'dz 
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where the distance from the Galactic center to the sun is taken to be 8 kpc and z denotes 
the height above Galactic disc in unit of kpc. This factor becomes important as the stars 
evolve to be luminous along the giant branch. For the dwarfs with low luminosity, on the 
other hand, the spatial variations in the density of EMP stars are negligible and we simply 
have Ns{L) oc L^/^. 

The size of surface convection, Mgcz, differs by a factor of ~ 100 (see eq. [16]) between 
dwarfs and giants, and hence, the upper limit, Am{Mi, [C/H], Mgcz), to the binary separation 
of carbon enhancement for dwarfs and subgiants is 10 times larger than that for giants. 
Thus, the number ratio of CEMP-s giants to CEMP-s dwarfs and subgiants will be given by 
V's.giant/l^s.dwarf + V's.subgiant) , each of which is defined as; 
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By adopting the limiting magnitude of 15 mag, we find the resultant proportions of RGB+HB 

stars, subgiants, and dwarfs at 63 %, 10 %, and 27 %, respectively. In our samples, there are 
13 giants, 4 subgiants, and 4 dwarfs (62 %, 19 %, and 19 %, respectively) among CEMP- 
s stars if we define dwarfs as stars with the surface gravity, log (7 > 4.2, and subgiants 
as 4.2 > log (7 > 3.5, and giants (including HE stars) as log (7 < 3.5, respectively. The 
above estimate based on the binary scenario gives a reasonable agreement with the current 
observations. 



In the preceding section, we have demonstrated that the statistical features of CEMP-s 
St subclass are understandable in terms of our binary scenario. In this section, we 

take up the problems of the observed frequencies of CEMP-s and CEMP-nos stars among 
EMP stars and of the total number of EMP stars in an attempt to inquire into the initial 
mass function (IMF) of EMP stars and to give an insight into their nature. 

From our scenario, we can estimate the proportions of CEMP-s and CEMP-nos stars 
to EMP binaries. The number, t/'bmary, of all the EMP stars born in binary systems in 
a flux-limited sample may be evaluated in the same way as we have derived the CEMP-s 
population in the preceding section. We may divide the EMP binaries into three categories; 
in addition to 

(1) white-dwarf binary systems, so far discussed, in which the primary stars in the mass 
range of 0.8 Mq < Mi < M^p have already become white dwarfs and the secondary stars 
still undergo nuclear fusion, we define two more of 

(2) low- mass binary systems in which both of the members are in the mass range of < 0.8 Mq 
and still alive as nuclear burning stars, and 

(3) supernova binary systems in which the primary stars in the mass range of Mi > M^p 
have already exploded as supernovae and the secondary stars still remains on nuclear burning 
stages. 

Then the total EMP stars currently observed are given by the sum of these three categories 
as; 



7. Initial Mass Function and Binarity of EMP Stars 



■^binary = V'WDB + '0LMB + V'SN, 



(32) 



where 




(33) 
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Here 74inin_wD(^i, -^2), ^mm,Ms(^i; ^2) and Ajnm,SN{Mi, M2) denotes the minimum sepa- 
rations that can contain two stars of AGB and main sequence, both of main sequence stars, 
and of supernova progenitor and main sequence star, respectively: A^^t denotes the upper 
bound of binary separation, and is taken to be consistent with P — 10^° days given by 
Duquennoy & Mayor (1991). 

The evaluation of these numbers requires the initial mass function of EMP stars in the 
whole mass range, and we may well assume a lognormal form of IMF with the medium mass, 
Mmd, and the variance. Am, as parameters, i.e., 

eWo.iexpf- "°«"'-.y"^n . (36) 



m 



2x 



For the Galactic spheroid, the IMF is shown to be well approximated by a lognormal form 
with the parameters of M^^ = 0.22 Mq and A^ = 0.33 for low mass stars (m < 1 Mq) and 
Salpeter form for more massive stars (see e.g., review by Chabrier 2003). EMP stars may not 
necessarily follow the same IMF as the Galactic spheroid components of larger metallicity. 
Rather we here treat these two parameters as free and discuss the constraints imposed on 
them from the comparisons with the current observations. As for the binary mass ratio and 
period, we assume the same distributions as adopted in the preceding section. 



7.1. Frequency of CEMP-s Stars 

The population of CEMP-s stars, ■0cemp-s, predicted from our binary scenario, is given 
by the sum of CEMP-s dwarfs, subgiants and giants (including HB stars) in eqs. (29)-(31); 

V'CEMP-s = V's.dwarf + 'V's.subgiant + "^s.giant- (37) 

We may integrate eqs. (37) and (32) by use of IMF in eq. (36). Figure 10 plots the proportion 
of CEMP-s stars to EMP binaries, ■0CEMP-s/'0binary, as a function of medium mass, Mmd, 
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with a fixed dispersion of Am = 0.33, assuming simple lognormal IMFs. If the same initial 
mass function as the Galactic Spheroid (peaked at 0.22 Mq) is applied, we find 

V'CEMP-s/^WDB - 0.58 : 
^CEMP-s/^binary - 0.14. 

The former ratio indicates that more than a half of white-dwarf binaries have produced 
CEMP-s stars, which is attributable to the very large coverage in the parameter space of 
separations and mass ratios, as stated above. Nevertheless, their proportion to the total EMP 
stars decreases by a factors of ~ 4.1 owing mainly to the contribution from the low-mass 
binaries. If we take into account the contribution of the stars born as single, this fraction 
is too small when compared with the observed very high frequencies of CEMP-s stars (see 
below). This may be indicative of more massive IMF for EMP stars. As Mmd increases, the 
frequency of low-mass binaries diminishes and the CEMP-s proportion increases to attain 
at a maximum of '0cEMP-s/V'binary — 0.37 near — -Mne-FDDM- For a still larger Mmd, it 
turns to decrease as the contributions augment from white-dwarf binaries with the massive 
primaries of M > Mhc-fddm and then from supernova binaries of M > Mup. In this figure, 
we also plot the CEMP-s proportion if we assume that the same number of EMP stars are 
born as single. The contribution of single-born EMP stars may be significant for the IMF of 
small Mmd, but decreases rapidly to be negligible for the IMF of large Mmd > Mne-FDDM- 

From the HK and HES observations, it is reported that the proportion of CEMP stars 
of [C/Fe] > 1 amounts to 25% for in the metallicity range of [Fe/H] < —2.5 (Beers 1999; 
Rossi et al. 1999; Christheb 2003, see also Lucatello et al. 2005a,b). They include the 
contribution from CEMP-nos stars, which accounts for ~ 25% of CEMP stars (Ryan et 
al. 2005; Aoki et al. 2006). If allowance is made for the contribution of CEMP-nos, then, 
the CEMP-s proportion may be smaller, while larger if we include such CEMP-s stars as 
having smaller carbon enrichment of 1 > [C/Fe] > 0.5. Recently, Cohen et al. (2005) find 
the frequency of CEMP stars at 14.4 ± 4% for [Fe/H] < -2 and [C/Fe] > 1.0 in HES 
sample, by correcting errors in abundance analyses of CEMP stars by means of the follow- 
up spectroscopy. This gives the CEMP-s fraction at ~ 10% with the possible contribution 
from CEMP-nos subtracted. However, their definition of CEMP stars is difTerent from ours 
and may underestimate the CEMP-s fraction since the efficiency of producing the carbon 
enrichment greatly decreases for larger metallicity of [Fe/H] > —2.5, as discussed above. 

Accordingly, we may well adopt the observed proportion of CEMP-s stars among the 
EMP stars at 10 ~ 25%. As seen from Fig. 10, there are two possible solutions with the 
separate ranges of Mmd that predict the proportion of CEMP-s stars compatible with the 
observations; one is the low-mass IMF with Mmd — 0.6 ~ 2.5 Mq and the other the high- 
mass IMF with Mmd — 7.5 ~ 14 Mq. The two solutions are separated by the overproduction 



-29- 



of CEMP-s stars from the white dwarf binaries because of large coverage in the parameter 
space, as seen above from Fig. 8. Although the both IMFs produce the similar proportions 
of CEMP-s stars, they predict quite different constituents for other EMP stars. In the case 
of the low-mass IMF, EMP stars other than CEMP stars are constituted mainly of low-mass 
stars born as single and of the members of low-mass binaries. In the case of the high- 
mass IMF, the extant EMP stars were most born as the low-mass member of binaries and 
more than a half of them stem from the binaries of massive primary, which have undergone 
supernova explosions. 

Our result differs from that of Lucatcllo ct al. (2005b), who find only the low-mass 
solution with Mmd — 0.79 M©, which arises from the different assumption on the distribution 
of mass fraction n{q); we adopt a flat distribution of mass fraction while they assume that 
the primary and secondary stars both obey the same initial mass function of single stars. 
As stated above, either of these two assumptions is not ruled out according to the current 
researches, and yet, they entail distinct conscqTicnces on the mass spectra of secondary stars. 
We have a nearly flat mass function for CEMP and also for CH stars. On the other hand, if 
we apply their assumption to CH stars, it entails the predominance of low mass stars down 
to ~ 0.2 Mq (Chabrier 2003). From the observation of CH stars, however, the contrary is 
reported that subgiant CH stars are found only in a narrow spectral-type range around GO 
but not among dwarfs of late G and K- types (Luck & Bond 1982). This lack of low-mass 
main-sequence CH stars favors our assumption of a flat distribution of mass ratio rather than 
the same single-star mass-function of primary and secondary stars. Luck & Bond (1991) have 
argued the thick accreted layers as the cause, but it seems difficult to prevent the formation 
of low mass CH stars through the wind accretion, as seen in § 5. These assumptions are 
distinguished also in the production of CEMP-nos, as discussed below. 



7.2. Relative Frequencies of CEMP-s and CEMP-nos Stars 

In our scenario, both CEMP-s and CEMP-nos arise from the white dwarf binaries, for 
which the existence and absence of s-process element enhancement are separated by the mass 
of primary components at Mi = Mhc-fddm- Then, the populations of CEMP-nos stars can 
be estimated in the similar way to the populations of CEMP-s stars as; 

'0CEMP-nos = 'ip nos, dwarf ~l~ V^nos, giant ~l~ V^nos, subgiant ; 

(38) 

where the each term in the right-hand side member can be computed by integrating the 
corresponding equations of (29)-(31) for the mass range of Mi between Mhc-fddm and M^p 
and by replacing the lower bound of separation by ^tdu- Because of difference in the 
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primary mass, the ratio between ■0cemp-s and ■0CEMP-nos depends on the assumed initial 
mass function. If we assume the Salpeter IMF for M > 0.8 Mq, it reduces to 

'0CEMP-nos/^CEMP-s - 1/50, 

which predicts a ncghgible fraction of CEMP-nos stars as compared with CEMP-s stars. 
The distribution, n{q), of mass fraction is unhkely to differ by an order of magnitude in the 
mass range of concern here and the distribution, /(-P), of orbital period is also expected not 
to vary greatly with the primary mass. Accordingly, larger population of CEMP-nos stars 
is possible only for the initial mass function shifting toward higher mass. 

In Fig. 10, we plot the proportion of CEMP-nos stars to the EMP binaries and the 
ratio between the CEMP-nos and CEMP-s stars, as a function of Mmd- As Mmd increases 
and the IMF shifts to be more massive, the CEMP-nos proportion first increases, and hits a 
maximum of V'CEMP-nos/^binary — 0.17 near ~ Mup. The maximum fraction is smaller 
than that obtained for CEMP-s stars above because of larger mass ratio between the primary 
and secondary stars. For a still larger Mmd, the CEMP-nos fraction turns to decrease and the 
supernova binaries outnumber the other type binaries. The ratio, V'CEMP-nos/V'CEMP-s) is a 
monotone increase function of Mmd! it starts from a very small values for the Salpeter IMF, 
the increase is accelerated for Mmd > ^Hc-fddm because of the decrease in the CEMP-s 
proportion, and eventually, the CEMP-nos proportion exceeds the CEMP-s proportion for 

M^d > 11 Mq. 

The ratio between CEMP-nos and CEMP-s stars is reported to be ~ 1/3 if limited to 
[C/Fe] > 1 (Aoki et al. 2006 in preparation; Tsangarides et al. 2004). In our compilation of 
[C/Fe] > 0.5, the ratio between CEMP-nos to CEMP-s stars is 13/21 = 0.62, with two stars 
rich in r-process elements included. In the above estimates, the upper mass bound, M^p, 
of primary stars for the binaries producing CEMP-nos stars is set equal to the lower mass- 
limit at which carbon ignite under the non-degenerate condition. Because of uncertainties 
in modeling of hot bottom burning (e.g., due to dependence on the mixing length), however, 
the boundary may not be clearly dehneated from those yielding EMP stars that have low 
or moderate carbon enhancement but show large nitrogen enrichment. In particular, most 
of the "mixed" stars, discussed by Spite et al. (2005) are left out from our compilation since 
we define CEMP stars as [C/Fe] > 0.5. These nitrogen-rich stars could be considered as 
the counterpart of CEMP-nos stars with more massive primary, and hence, processed deeply 
by the hot bottom burning. If "mixed" stars are included, the number of CEMP-nos stars 
amounts to 27 and is comparable with, or even exceeds, the number of CEMP-s stars in 
our list. Accordingly the current observations suggest the ratio of CENP-nos to CEMP-s 
stars in the range of V'CEMP-nos/V'CEMP-s ^1/3 and we may well take an upper bound at 
~ 1. The ratios in this range agree well with those predicted from the high-mass IMF with 
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Mmd ~ 10 M0, one of the two solutions derived above from the CEMP-s proportion. The 
low-mass IMF is excluded on the basis of current observations since it can predict too small 
proportion of CEMP-nos stars ('0CEMP-nos/'0CEMP-s < 0.05). 

7.3. Nature of BMP Stars in Galactic Halo 

We have shown that the IMF for EMP binaries can be constrained from the observed 
characteristics of CEMP-s and CEMP-nos stars. The constraints thus derived are summa- 
rized in Figure 11 on the diagram of the medium mass Mmd and the standard deviation Aa/ , 
the parameters of initial mass function in the lognormal form of eq. (36). Solid lines denote 
the loci of the IMFs that yield CEMP-s stars at constant fractions, and broken lines the loci 
of IMFs that produce CEMP-s and CEMP-nos at constant ratios. Since the production rate 
of CEMP-s stars from the white dwarf binaries is very high (~ 60%), the parameter space 
near M^^d — ^Hb-fdom is excluded by the overproduction of CEMP-s as compared with the 
observation for the dispersion of Am ^ 0.4. For smaller dispersion, therefore, the allowed 
parameter space for CEMP-s stars is separated into two parameter spaces with low and high 
mass regimes, converging toward ~ 1 and 6Mq respectively. As Am increases, the 

wider range of comes to be compatible with the observation of CEMP-s stars, and two 
regimes are connected for such broad IMF as Am > 0.4. On the other hand, the formation 
of CEMP-nos stars at the observed ratios to CEMP-s stars restricts the IMFs to those with 
the large Mmd (> Afne-FDDM); the allowed range of Mmd increases as the IMF grows broader 
and gives birth to larger fraction of low mass stars. 

In conclusion, only compatible with the observations both of CEMP-s and CENP-nos 
stars are the IMFs with Mmd greater than ~ 6 Mq and increasing with Am, as shaded in 
Fig. 11. One of the important consequences is that the currently observed EMP stars have 
to stem exclusively from the binary systems; Figure 12 exemplifies the mass distributions 
of primary and secondary stars for a typical initial mass function derived here. We see 
that the secondary components well extend into low mass regime below 0.8 Mq but the 
primary components have a negligible fraction of low mass stars. We may not expect a 
significant fraction of low mass stars born as single unless the initial mass function for single 
stars may greatly differs from that of primary stars. Another remarkable consequence is 
that a significant fraction of currently observed EMP stars were formed as the member of 
supernova binaries of Mi > M^p. For example, for M^^ in the range obtained in Fig. 10, the 
proportion of supernova binaries accounts for 40 ~ 60% of the currently observed EMP stars. 
The proportion of low-mass binaries is less than 0.1 % and the contribution to EMP stars 
is less then 1 %. Note that our conclusion is dependent on the assumption of n(q), which is 
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presently very uncertain. In order to explain the origins of not only CEMP-s stars but also 
CEMP-nos stars, however, the IMFs of EMP stars have to be weighted in the intermediate 
and more larger mass range. Otherwise, one has to seek other formation mechanism(s) for 
all the CEMP-nos stars, which are utterly unknown to the current theory of stellar evolution 
and/or nucleosynthesis. 

The present results further provide a way of probing into the stellar populations that 
have left EMP stars now constituting the Galactic halo. Wc may estimate the total stellar 
mass necessary to explain the number of currently observed EMP stars. From the above 
derived IMF, wc expect one low-mass star of M < 0.8 Mq out of EMP binaries of number 
A^b = -^md/0.8 Mq, and hence, of total stellar mass (3/2)A^b^md on average if we assume the 
same flat distribution of mass ratio [n{q) = 1] as above. On the other hand, from the recent 
large scale survey, the number of observed EMP stars in the Galactic halo is estimated at 
~ 670 steradian"^ for [Fe/H] < —2.5 and with the hmiting magnitude oi B < 17.5 (Beers & 
Ghristlieb 2005, in deriving this number we take the ratio between EMP stars of [Fe/H] < —2 
and [Fc/H] < —3 from their Table 3). With this limiting magnitude, giants can be reached 
to distance up to ~ 100 kpc, and hence, within the whole stellar halo, while dwarfs can be 
reached only in neighborhood of ~ 3 kpc. Giants occupy about a half among the observed 
EMP stars, as discussed in § 4. By taking into account the mass range of stars now on the 
giant branch, AMq = 0.01 M© and the flat mass function of EMP stars, we may estimate 
total number of EMP stars in the Galactic halo, 

Nemp ^ 670 X 0.5 X 47r X 0.8 Mq/AMq ~ 3.4 x 10^ (39) 

Thus, the total mass, M[Fe/H]<-2.5, of stars in the mother stellar populations of [Fe/H] < —2.5 
that have produced these low-mass EMP stars that survive to date in the Galactic halo, 
amounts to be 

-^[Fe/H]<-2.5 ~ 6 X 10^ Mo (7Vemp/3 • 10^) {M^JW M^f . (40) 

The loci of constant mass of M[Fe/H]<-2.5 (= 10^, 10^, and 10^ Mq), obtained numerically 
from the equations given in the preceding subsections, are plotted in Fig. 11 (dash dotted 
lines). The total mass of stars in the mother populations increases with Mmd, and most of 
them have exploded as supernovae. Accordingly, the metal production by these erstwhile 
supernovae may impose an upper bound on the total stellar mass that has been involved in 
the mother populations, and hence, an upper bound on the medium mass. If we take the 
averaged iron yield to be ~ 0.01 of the initial stellar mass, then, the mother populations 
may have increased the iron abundance in our Galaxy up to 



[Fe/H] ~ -2 + log (M[pe/H]<-2.5/10' Mq) 



(41) 



-33- 



on an averaged basis all over the Galaxy of total (baryon) mass M ~ 10^^ Mq. With the 
derived high-mass IMF, the total mass of stars of 10^ Mq are sufficient to promote the chem- 
ical evolution of the Galaxy to the metallicity of population II. As a consequence, the IMFs 
with significantly exceeding ~ 10 may be excluded by the metal overproduction. 

In summary, the EMP stars currently observed originate from a small fraction (~ 10% 
of binary systems in number) of stellar populations of total mass ~ 10^ Mq that have once 
constitute, or merged into the Galactic halo. A significant portion of them become CEMP 
stars, and other significant portion (40 ~ 60%) have been exposed to supernova explosion 
of their companions. In the supernova binaries, the secondary stars arc likely to be set 
unbound after the supernova explosions as a result of sudden reduction of the primary mass. 
The secondary stars may possibly interact with the envelope matter lost through the wind 
before being released, and also with SN ejecta from the primary stars. 

We conclude this section with two comments on the implications of these supernova 
binaries to the observed characteristics of EMP stars. Because of large wind velocity and 
the expansion velocity of supernovae, only a small fraction of ejecta can be accreted by the 
secondary stars, and yet, it may influence the surface characteristics for the elements of such 
small abundances as the r-process elements. The r-process is most poorly understood among 
the stellar nucleosynthesis mechanisms, and yet, it may be agucd from the solar abundances 
and the supernova rates that the amount of r-process elements, ejected per one event, is of 
an order of M,_p ~ 10"^ Mq on average (Mathews & Cowan 1990; Woosely et al. 1994). 
Then, simply assuming the geometrical cross sections for the accretion, we may expect the 
surface enrichment of r-process elements of secondary stars as large as 

[r/Fe] ~ 1.3 + log(M,_p/10-^ Mq) - 21og^(AU) 

- log(M,ez, g/0.35)Mq (42) 

for the giants of [Fe/H] ~ —3. EMP stars are known to display variations of r-process 
element abundances with a large range by a factors of ~ 1000 of —1 < [Eu/Fe] < 2 (Honda 
et al. 2004). The above estimate can be compatible with the observed enrichment, with the 
largest one from the systems of the smallest separations. In particular, our binary scenario 
gives a straightforward explanation to the observed large variations in term of the difference 
in the binary separation. This new channel of surface pollution is worth future investigations 
with the interactions between the matter ejected by supernova explosion and the secondary 
star taken into account. It may also happen that the secondary stars accrete the envelope 
mass ejected by wind from primary stars before the supernova, and also, be polluted through 
the accretion in gas shell of supernova remnants after the explosion. 

Among the supernova binary systems, the primary stars of M^p < Mi < 11 Mq have 
been proposed to make supernova explosion, triggered by electron capture on ^°Ne and 
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oxygen burning in the electron-degenerate O + Ne core by Miyaji et al. (1980). On the 
other hand, Ritossa et al. (1996) show that the helium layer is dredged up by the surface 
convection during the carbon shell-burning, and then, carbon shell-burning is extinguished 
and these stars enter into the thermal pulsating (Super) AGB phase with the hydrogen 
and helium double shell burnings. If this is the case, these stars end in O + Nc white 
dwarfs, ejecting the envelope through wind mass loss, just as AGB stars with C -|- O core 
(Ritossa et al. 1999; Gil- Pons et al. 2005), although lower efficiency of mass loss may tend to 
narrow the mass range for stars of lower metallicity. Consequently, some of the binaries with 
the primary stars in this mass range may produce EMP stars of small carbon abundance 
but largely enriched with nitrogen, similar to "mixed" stars, because of dilution in larger 
envelope mass and of deeper processing by the hot bottom burning on larger core mass at 
the onset of TP-SAGB phase, as in the upper mass end of white dwarf binary systems of 
M < Mup. In addition, because of the very short lifetimes of these massive primaries (< 10® 
yr), the secondary stars may suffer from the surface pollution by later accreting interstellar 
gas, enriched with metals, if their parent clouds persist sufficiently long to be polluted by 
supernova ejecta of subsequent generations (Suda et al. 2004). 



8. Conclusions and Discussion 

In this paper, we propose that the origins of carbon-enhanced metal-poor (CEMP) stars, 
currently observed in the Galactic halo, are explained in terms of the evolution of binary 
systems of extremely metal-poor (EMP) stars of [Fe/H] < —2.5. Our binary scenario is 
based on the evolutionary models for both EMP stars and binary mass transfer. We have 
examined it by confronting the consequences with the observations such as the relation of 
the carbon-enhancement to the orbital periods and the statistical properties of CEMP stars 
and also through the comparison with their counterparts of higher metallicity, CH stars. 
Finally, we demonstrate that a constraint can be imposed on the initial mass functions of 
EMP stars and discuss the implications to nature of EMP stars in the Galactic halo. 

Main results are summarized as follows: 

(I) We present an updated summary of cvohition of EMP stars. The primary stars can 
develop the surface carbon enhancement via two distinct mechanisms when evolving to AGB; 
He-flash driven deep mixing (He-FDDM) in the mass range of M < 3.5 Mq and third dredge- 
up (TDU) in the mass range of 1.5 Mq ^ M < Mup, where M^p denote the upper mass limit 
of stars ended as white dwarfs. Nitrogen enhancement results either from He-FDDM or from 
the hot bottom burning (HBB) in the envelope, and s-process elements are synthesized in 
the helium convective zone promoted by hydrogen mixing precedent to He-FDDM. 
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(II) The secondary stars accrete a part of carbon-rich envelope, ejected from the AGB 
primary stars, to be CEMP stars. CEMP-s stars with nitrogen and s-process enhancement 
originate from the EMP binary systems with the primary stars of mass between 0.8 < 
^1 ^ 3.5 Mq through He-FDDM. On the other hand, CEMP-nos stars without s-process 
enhancement stem from the systems with the primary of mass 3.5M0 < Mi < M^p. In the 
latter case, the abundance of nitrogen as well as of carbon vary with the mass of primary, and 
nitrogen-rich stars with mild or little carbon enhancement are associated with the massive 
primaries. 

(III) Our binary scenario is shown to give reasonable accounts to the observed characteristics 
and statistical features of CEMP stars. The large fraction of CEMP stars are explained by 
the broad parameter space of mass range and binary separation of progenitor binary systems 
as compared with CH stars. In particular, it is shown that CEMP binary systems have the 
orbital periods of much wider range than CH stars; all the detected orbital periods and 
the lack of detected variations in the radial velocity from about a half of CEMP stars are 
both consistent with the predictions from the binary scenario and from the time span of 
monitoring of binarity to date. 

(IV) Prom the observed frequencies of CEMP-s and CEMP-nos stars among EMP stars, we 
demonstrate that the initial mass function (IMF) of EMP stars has to be massive with the 
medium mass, Mmd ^ 6 Af© if approximated to a lognormal form. This also implies that the 
currently observed EMP stars in the Galactic halo were formed exclusively as the members 
of binary systems; low-mass EMP stars born as single stars account for a tiny fraction (less 
than ~ 1%) if the single stars were born at nearly equal numbers with the similar initial 
mass function to that of primary stars. Accordingly, in addition to CEMP stars now with 
white dwarf companions, a significant fraction of EMP stars (40 ~ 60%) used to have a 
massive primary star and have been exposed to supernova explosion before dismissed from 
the binary systems. We suggest, as a new channel of surface pollution of EMP stars, the 
accretion from stellar wind matter and SN ejecta of massive companions. 

(V) Prom the total number of EMP stars in the Galactic halo, obtained from the recent 
large-scale surveys, the total mass of stars in their mother stellar populations is estimated to 
be ~ 10*^ M0(Mmd/lO Mo)^. The metal production of SNe imposes another constraint and 
may rule out the very high mass IMF with the medium mass, Mmd, significantly exceeding 

10 Mq. 

The following picture of Galactic halo emerges from the present study. Galactic halo 
has once involved the stellar populations of metallicity [Fc/H] < —2.5 and of mean mass 
~ 10 M0, which contain the binary systems as many as ~ 10^(Mind/10 Mq) in total; about 
a half of them had the primary components of low and intermediate masses that have now 
evolved to white dwarfs, and another half have more massive primary components that have 
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ended as supernovae. The currently observed EMP stars are ~ 10% survivors of the low-mass 
members of these binaries. The mother stellar populations of EMP stars therefore constitute 
0.1% of baryon mass of our Galaxy (~ 10^^ Mq). SNe from the massive stars comprised in 
them suffice to raise the metallicity of whole Galaxy to ~ 0.01 solar on average, leading to 
the formation of population II objects, observed in the Galactic halo. 

Our results indicate that the transition from an IMF dominated by massive stars to an 
IMF overwhelmed by low- mass stars occurred after the metallicity is raised above [Fe/H] ~ 
—2.5. Our interpretation on the origin of CEMP-nos stars, i.e., TDU and hot bottom 
burning, should be common to the stars of more metal rich generations. Therefore, the lack 
of their correspondences in Pop II stars, particularly nitrogen-rich "mixed" stars, may be 
attributed to the difference in the initial mass function, as discussed in §7. The metallicity 
at the transition suggested from the present work seems significantly larger than claimed 
from the studies on the dynamical and thermal evolution of gas clouds under the metal 
deficient circumstance. The recent studies tend to argue that the metal and dust cooling can 
supersede the cooling by hydrogen molecules even for the metallicity as small as [Fe/H] ~ —5, 
which is argued to reduce the Jeans mass and produce the fragments below a solar mass 
(e.g., see Omukai 2000; Omukai et al. 2005). It is also shown that the sub-solar Jeans 
masses and fragments can be formed even from gas completely devoid of metals, once gas 
is heated to be ionized by a shock (Uehara & Inutsuka 2000) either due to the collapse of 
massive primordial objects with total (baryon plus dark matter) mass M > 10^ Mq or due 
to supernova explosion of the first generation stars (Machida et al. 2005). In these studies 
the Jeans mass is directly connected to the masses of formed stars. It is true however that 
the process of fragmentation is still poorly understood and a proper understandings of star 
formation is yet to be established. It is worth noticing here that the globular clusters, which 
embrace a host of low-mass stars, exist only for the metallicity of [Fe/H] > —2.5 in our 
Galaxy. 

One of the critical assumptions in the present study is on the distribution of mass ratio 
n{q) of binary systems. It is still an open question either observationally or theoretically, 
as discussed in §6.1. We postulate a fiat mass function of secondary components extending 
well below a tenth of the mean mass of primary components, which enables us to explain the 
origins not only of CEMP-s stars but also of CEPM-nos stars within the current standard 
framework of the theory of stellar evolution. If both the binary components are assumed to 
have the same distribution functions, as done by Lucatello et al. (2005b), we should seek 
for the other origin of CEMP-nos stars presumably elsewhere outside the current theory. 
In addition, we have pointed out that a flat distribution of mass ratio finds support in the 
low- mass cut off reported for CH stars (Luck & Bond 1991). 
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Another critical assumption in our argument is on the s-process nucleosynthesis in EMP 
star. In the present paper, we work out the stellar evolution and nucleosynthesis within the 
standard frameworks, which takes into account the thermal convection and chemical diffusion 
as the mechanisms of material mixing in the stellar interior with a negligible contribution 
of convective overshooting. It is true that there exist the observations of surface abundance 
anomalies, not explicable within these frameworks. One of the relevant issues is the formation 
of ^^C pocket, which are proposed to work as neutron sources in the metal-rich stars, but 
can not be realized within the current standard framework of stellar evolution. Our scenario 
postulates that it is inefficient and the s-process nucleosynthesis via radiative ^^C burning 
is not effective in the stars of metallicity below [Fe/H] < —2.5, at least, in the massive 
AGB stars of M > 3.5 Mq. In actuality, the necessity of ^^C pocket is claimed from the 
comparisons with the observed distribution of s-process elements in stars of more metal 
rich populations, and yet, it is treated as a free parameter. The present results may give 
an insight into the modeling of ^^C pocket formation if our interpretation of CEMP-nos is 
correct. 

As for s-process nucleosynthesis in EMP stars, we may point out that the relevance of 
our high-mass IMFs to the synthesis, in particular, of the light s-process elements. Most of 
CEMP-nos stars in our samples show larger abundances of strontium than barium ([Sr/Ba] > 
0) , distinct from CEMP-s stars which show larger enrichment of the main and heavy s-process 
elements relative to the light s-elements. This may be explicable in terms of the difference 
in the mass of primary stars, and hence, in the core mass when the primary stars start 
the TP-AGE evolution. For the primary stars of CEMP-nos stars, the core mass can be 
sufficiently large and the temperature in the helium flash convection high enough to burn 
^^Ne via ^^Ne(a, n)^^Mg, and to promote the s-process nucleosynthesis (Iben 1975b). This 
neutron capture process is likely to produce mostly light s-process elements and to yield such 
distribution with the ratio of [Sr/Ba] > 0, as shown by Truran & Iben (1977), for the very 
source reaction produces neutron poison, and hence, the available neutron per seed nuclei 
are restricted to be rather small. Recently Aoki et al. (2005) report a general tendency 
of the excess of light over main s-process elements increasing toward lower metallicity of 
[Fe/H] < —2.9 up to [Sr/Ba] ~ 1.5. The hehum shell flashes in massive AGB stars can be 
one of the candidates for the site of hght s-process element synthesis. 

As for some CEMP stars, there have been proposed other origins than the binary evo- 
lution; Umeda & Nomoto (2003) construct a peculiar supernova model with a carbon-rich 
ejecta and propose the second generation of stars formed of gas mixed with the ejecta as a 
possible scenario for CEMP stars. In particular, for stars below [Fe/H] < —5, recently dis- 
covered, the formation of low mass stars is argued by assuming the high carbon abundance 
[C/H] ~ —1.3 as observed from the stars (Bromm & Loeb 2003). Since the amount of carbon 
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yield from type II SNe is of ~ 0.2 M©, however, it is open to question whether the second 
generation stars were formed of carbon abundance as large as [C/H] = —1.2 ~ —1.3 with 
such a small amount of carbon ejected. This seems at an apparent variance to the fact that 
most of EMP stars have iron abundances of [Fe/H] < —2.5 although they are thought to be 
formed with iron ejecta of mass 0.1 ~ IMq in the similar way; Further with the smallest 
energy of explosion < 10^^ erg assigned to the carbon-rich SNe, the star formation triggered 
by supernova (Tsujimoto et al. 1999) is unlikely to work since the shock driven by SN dies 
away before the fragmentation occurs in the swept shell (Machida et al. 2005). Even if the 
stars of the second generation are not formed, however, our binary scenario may propose an 
alternative channels of the surface pollution of the low-mass members of the first generation 
by accreting the wind matter and supernova ejecta. This is worth consideration in future 
works. 

Finally, we comment on the surface pollution by accreting interstellar gas in the parent 
clouds where these EMP stars were born, discussed by Shigcyama et al. (2003) and by Suda 
et al. (2004). The latter argue that the surface of EMP dwarfs may be polluted up to 
[Fe/H] > —3 if the parent clouds survive sufficiently long (> 10^ yr) and enrich their gase 
with metals up to mctallicity of [Fe/H] ~ —2 and higher. If the lifetimes of primary stars 
are shorter than < 10^ yr as in the cases of white dwarf binaries at the massive end and of 
supernova binaries, EMP stars may have suffered from the accretion of interstellar gas after 
the mass transfer, and hence, can disguise their surface abundance with that of interstellar 
gas. When they evolve to giants, however, the surface abundance is diluted by a factor of 
~ 100 with the internal matter, accumulated by the mass transfer from the primary stars. 
Accordingly, some of CEMP-nos stars and "mixed" stars, and some from the supernova 
binaries, may have quite different appearance while they are dwarfs and after they evolve to 
giant branch. 

The surface pollution by accreting interstellar metal-rich gas has been addressed out 
in relation to the two most iron-deficient, carbon-enhanced stars of [Fe/H] < —5 recent 
discovered by Christhcb et al. (2002) and by Frcbel ct al. (2005). Suda et al. (2004) have 
suggest the possibility that these stars are Pop 111 stars which have been polluted by accreting 
iron-rich interstellar gas, and also, have become CEMP stars through the mass transfer from 
the erstwhile AGB companion in the close binaries. If they are really Pop III survivors, 
we may expect that some Pop III stars without carbon enhancement exist, too. These 
two stars are thought to undergo He-FDDM, and hence, if the Pop III stars were formed 
under the same IMF as obtained above, then, we may expect 8 ~ 20 Pop III stars mostly 
without carbon enhancement (3 ~ 12 single stars from the supernova binaries). Since their 
iron group elements are due to pollution with mctal-rich interstellar gas, these stars have 
to display the surface metallicity of [Fe/H] ~ — 3 ~ — 4 (with the gravitational settling is 
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taking into account) while they are dwarfs and suffer from the dilution of surface pollution 

as they evolve to giants and deepen the surface convection. The surface abundances of these 
pollutants should be similar to those observed from HE0107-5240 and HE1327-2326 for the 
elements accreted from the interstellar gas since for metallicity of [Fc/H] ~ —2 and higher, 
the variations in the abundances tend to be small. We may anticipate detection of such stars 
even current compilations of EMP stars with careful investigations. 

We are very much benefited from discussion with Dr. Icko, Iben, Jr. and Dr. Toshitaka 
Kajino. This work has been partially supported by Grant-in-Aid for Scientific Research 
(15204010, 16540213, 18104003), from the Japanese Society for the Promotion of Science. 

A. An attempt to Interpret the abundance anomalies from EMP stars 

A.l. CEMP-nos stars with Oxygen enhancement 

The hot bottom burning may burn not only carbon but also oxygen into nitrogen in 
metal-poor stars since the temperature in the bottom of surface convection increases with 
decreasing the metallicity (e.g., Ventura et al. 2001). If this is the case, the abundance of oxy- 
gen now observed should be inherent in the secondary stars, mixed and diluted with matter 
transferred from the primary stars. The carbon abundances larger than the equilibrium ratio 
to nitrogen may also be attributed to the pristine matter of secondary stars. All CEMP-nos 
stars listed in our sample are giants and the accreted matter has been mixed in the deep 
convection. As a corollary, before developing the deep convection, these CEMP-nos stars 
should have a different appearance as dwarfs that only nitrogen is enhanced; in addition, the 
surface abundance of dwarfs may also be subject to the pollution by accreting interstellar 
matter after the mass transfer because of relatively short lifetime of primary stars, as dis- 
cussed by Suda et al. (2004). It is to be noted that the nitrogen enrichment in CEMP-nos 
stars is different from that observed from some giants in the globular clusters, attendant with 
the depletion of oxygen, and hence, the deep mixing mechanisms along the giant branch as 
proposed for the latter (e.g., see Suda & Fujimoto 2006) will not be applicable. 

A. 2. Nitrogen-rich EMP stars and the Relation to CEMP-nos Stars 

Some EMP stars show large nitrogen enhancement but a httle or no carbon enhance- 
ment. A well-known example is CD— 38°245, which were thought as the most metal-poor 
star until the discovery of HE0107-5240 ([Fe/H] = -4.5; Bessel & Norris 1984, but revised 
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later to [Fe/H] = -3.98 by Norris et al. 2002), with the abundances of [C/Fe] = - 0.3, 
[N/Fe] = 1.7, [0/Fe] = 1.3, and [Sr/Fe] ~ [Ba/Fe] = -0.5 (Bessel & Norris 1987, cf. 
Spite et al.2005). For comparison, we plot this star in Fig. 7 with two more similar stars, 
CS22878-101 ([Fe/H] = -3.25) and CS22952-015 ([Fe/H] = -3.43). These stars may arise 
from the same mechanism as CEMP-nos stars with nitrogen enrichment but may belong 
to the binaries with more massive primary and of wider separation, yielding greater pro- 
cessing by hot bottom burning and smaller enhancement of [C + N/Fe]. Spite et al. (2005) 
observe EMP giants with weak or no carbon enhancement and argue that they are grouped 
into "mixed" stars with [N/C] ~ 1, including above three stars, and "unmixed" stars of 
[N/Fe] < 0.5, both with similar [C + N/Fe]. They compile 17 "mixed" stars, among which 4 
stars have [C + N/Fe] > 0.5, and our list of CEMP-nos stars shares one star, CS22949-037, 
of [C/Fe] > 0.5. For the "mixed" stars, the problem is also reduced to identify the site(s) 
for conversion of carbon into nitrogen, and they suggest several possible sites, including the 
hot bottom burning in massive AGB stars and the hydrogen burning in the envelope of 
very massive stars. In the latter case, however, we have to seek a way to retain such large 
abundance ratios of nitrogen to carbon and to oxygen as observed from the "mixed" stars if 
these stars were formed from matter polluted with matter ejected by supernova explosions 
of the massive stars since the nitrogen-rich matter from their envelope suffer from dilution 
due to the mixing with carbon- and/or oxygen-rich cjccta from the inner parts as well as 
with the interstellar matter (but see Norris et al. 2002). Instead, the origin of "mixed" stars 
may well be understood in terms of the binary evolution of case IV' but with the primary of 
large masses and of large separations. 



A.3. Abundance Features of CS22892-052 and Relation with Orbital Period 

CS22892-052 is of metalhcity [C/H] ~ —2 and observed to show the orbital period of 
P = 127 days. In order to explain the origin in terms of binary scenario, this star seems to 
require the small carbon abundance [C/H]p < — 1 in the wind, as seen from the comparison 
with the plots in Fig. 9. The binarity of this system is subject to suspicion, however, since 
it displays radial velocity variations of only a very low amplitude (~ 1 km s^^; Preston & 
Sneden 2001). If this star really belong to a binary of P = 127 days, the low amplitude 
implies either a significantly small inclination angle of the orbital plane and/or a very small 
mass ratio. In the latter case, the companion star is likely to be a brown dwarf rather than a 
white dwarf, and we have to seek the origin of carbon and other heavy elements of this star in 
the other site such as the preceding supernovae. More specifically, Tsujimoto & Shigeyama 
(2001) have argued that this star with enhanced abundances of r-process elements inherited 
its heavy elements from the supernova explosion of a ~ 20 star. On the other hand. 
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this star shows the enhancement and anomalous relative abundances of CNO abundances, 
as seen from Fig. 7, and this makes the mass transfer in a binary more likely as long as CN 
enrichment is concerned. 
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Table 1. Observational Sample of CEMP Star from Literature 



Object Period e log g [Fe/H] [C/Fe] [N/Fe] [0/Fe] [Sr/Fe] [Ba/Fe] [Pb/Fe] Ref. 



stars with radial velocity variation 



HE0024-2523 


3.14 





4.3 


-2.72 


2.6 


2.1 


0.40 


0.34 


1.46 


3.3 


1 


G77-61 


245 





5.05 


-4.03 


2.6 


2.6 


/ 





<1 


/ 


2 




342 


n 


4 1 


-2 57 


2 30 


2 12 


1 48 


84 


2 17 


3 55 


3 4 5 


CS22948-027 


426.5 


0.02 


1.8 


-2.47 


2.43 


1.75 


/ 


0.90 


2.26 


2.72 


6,7,8 


CS22942-019 


2800 


0.1 


2.4 


-2.64 


2.0 


0.8 


/ 


1.7 


1.92 


<1.6 


8,9,10 


CS22957-027 


3125 


0.45 


2.4 


-3.11 


2.4 


1.6 


/ 


-0.56 


-1.23 


/ 


8,9,10 


CS29497-034 


4130 


0.02 


1.8 


-2.90 


2.63 


2.38 


/ 


1.00 


2.03 


2.95 


7 

1 


LP625-44 


>12yr 


/ 


2.5 


-2.72 


2.25 


0.95 


1.85 


1.32 


2.81 


2.55 


11,12,13 ^ 


CS22892-052 


127.8? 




1.6 


-3.03 


0.89 


0.71 


0.72 


0.44 


0.92 


1.2 


8,13,14,15,161 


CS30301-015 


binary ^ 


/ 


0.8 


-2.64 


1.6 


1.7 


/ 


0.3 


1.45 


1.7 


9,10 


CS29526-110 


binary 


/ 


3.2 


-2.38 


2.2 


1.4 


/ 


0.88 


2.11 


3.3 


9,10 


HE2148-1247 


binary 


/ 


3.9 


-2.3 


1.91 


1.65 


/ 


0.76 


2.36 


3.12 


17 


CS22877-001 


binary? 


/ 


2.2 


-2.85 


1.0 


0.0 


/ 


-0.12 


-0.49 


/ 


6,18,19 


CS22183-015 


binary? 


/ 


2.5 


-2.85 


2.34 


/ 


/ 




2.09 


3.17 


19,20 


itars with s-process enhancement 




















HE1327-2326 


/ 


/ 


? 


-5.4 


4.1 


4.5 


<4.0 


1.1 


<1.4 


/ 


21 


HE0107-5240 


/ 


/ 


2.2 


-5.28 


3.70 


2.57 


2.3 


-0.52 


<0.82 


/ 


22,23 


HE0007-1832 


/ 


/ 


3.8 


-2.65 


2.55 


1.85 


/ 


0.14 


0.16 


<2.82 


24 


CS22947-187 


/ 


/ 


1.3 


-2.49 


1.03 


/ 


/ 


0.57 


1.18 


/ 


25 


CS31062-012 


/ 


/ 


4.5 


-2.55 


2.1 


1.2 


/ 


0.3 


1.98 


2.4 


9,10 


HD187216 


/ 


/ 


0.4 


-2.48 


1.3 


0.2 


/ 




2.5 


/ 


26 


CS31062-050 


/ 


/ 


3 


-2.31 


2 


1.2 


/ 


0.91 


2.30 


2.9 


9,10,27 



Table 1 — Continued 



Object 


Period 


e 


logg 


[Fe/H] 


[C/Fe] 


[N/Fe] 


[0/Fe] 


[Sr/Fe] 


[Ba/Fe] 


[Pb/Fe] 


Ref. 


CS22898-027 


/ 


/ 


3.7 


-2.26 


2.2 


0.9 


/ 


0.92 


2.23 


2.84 


9,10 


HD196944 


/ 


/ 


1.8 


-2.25 


1.32 


1.3 


/ 


0.84 


1.10 


1.9 


9,10 


stars with no s-process enhancement 


















CS22949-037 


/ 


/ 


1.5 


-3 9 


1 17 


2 57 


1 98 


33 


-0 58 


/ 


1 5 28 29 


CS29498-043 


/ 


/ 


0.6 


-3.75 


1.9 


2.3 


2.9 


-0.35 


-0.45 


/ 


10 


CS22891-200 


/ 


/ 


1 


-3 5 


53 


/ 
/ 


/ 
/ 


-1 33 

X . • ) • J 


-0 92 




23 


CS22897-008 


/ 


/ 


1.5 


-3.41 


0.56 


/ 


/ 


0.63 


-1.23 


/ 


15,24 


BS16929-005 


/ 


/ 


2.7 


-3.09 


0.92 


/ 


/ 


0.28 


-0.59 


/ 


14 


CS30325-094 


/ 


/ 


2.1 


-3 


0.5 


/ 


/ 


-2.35 


<-1.76 


/ 


30 


CS29516-041 


/ 


/ 


2.5 


-3 


0.5 


/ 


/ 


-1.95 


-2.00 


/ 


29,30 


CS30314-067 


/ 


/ 


0.7 


-2.85 


0.5 


1.2 


/ 


-0.37 


-0.57 


/ 


6 


BS16033-008 


/ 


/ 


2.7 


-2.8 


0.6 


/ 


/ 


-0.68 


-1.24 


/ 


30 


CS30325-028 


/ 


/ 


1.8 


-2.8 


0.6 


/ 


/ 


0.27 


-0.62 


/ 


30 


CS29502-092 


/ 


/ 


2.1 


-2.76 


1 


0.7 


/ 


-0.4 


-0.82 


/ 


6 



Note. — We select CEMP star as [Fe/H] < -2.5 and [C/Fe] > 0.5 from literature. Also included are N rich star with [Fe/H] < -2.2 and 
[C/Fe] > 0.5. 

''CS31062-012=LP706-7. 

'^Radial velocity variations are observed, but not enough to estimate period. 
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Table 2. Observational Sample of CH or Ba Star from Literature 



Object 


binarity 


Period (days) 


[Fe/H] 


[C/H] 


[C/Fe] 


[Ba/Fe] 


class 


Ref. 


HD77247 


y 


80.55 


/ 


/ 


/ 


/ 


Ba 


1 


HD46407 


y 


458.6 


-0.42 


0.11 


0.53 


1.39 


Ba 


1,2 


HD199939 


y 


584.9 


/ 


/ 


/ 


/ 


Ba 


1 


HD58368 


y 


672.7 


/ 


/ 


/ 


/ 


Ba 


1 


HD31487 


y 


1066.4 


/ 


/ 


/ 


/ 


Ba 


1 


HD223617 


y 


1301 


/ 


/ 


/ 


/ 


Ba 


1 


NGC2420-X 


y 


1402 


-0.7 


/ 


/ 


/ 


Ba 


1,3 


HD13611 


y 


1642.1 


/ 


/ 


/ 


/ 


Ba 


4,5 


HD101013 


y 


1710.9 


/ 


/ 


/ 


/ 


Ba 


1 


HD49641 


y 


1768 


/ 


/ 


/ 


/ 


Ba 


1 


HD 16458 


y 


2018 


/ 


/ 


/ 


/ 


Ba 


1 


HD204075 


y 


2422 


/ 


/ 


/ 


/ 


Ba 


1 


HD205011 


y 


2853 


/ 


/ 


/ 


/ 


Ba 


1 


HD178717 


y 


2866 





0.21 


0.21 


0.7 


Ba 


1,3,6 


HD131670 


y 


2948 


/ 


/ 


/ 


/ 


Ba 


1 


HD 196673 


y 


4000 


/ 


/ 


/ 


/ 


Ba 


1 


HD 199394 


y 


4390 


/ 


/ 


/ 


/ 


Ba 


1 


HD139195 


y 


5324 


/ 


/ 


/ 


/ 


Ba 


5,7 


HD202190 


y 


6489 


/ 


/ 


/ 


/ 


Ba 


5,7 


HD5825 


/ 


/ 


-0.7 


-0.85 


-0.15 




Ba 


8 


HD4084 


/ 


/ 


-0.7 


-0.58 


0.12 




Ba 


8 


HD15589 


/ 


/ 


-0.7 


-0.45 


0.25 




Ba 


8 


HD27271 


/ 


/ 


-0.5 


-0.35 


0.15 




Ba 


8 


HD44896 


/ 


/ 


-0.4 


0.36 


0.76 


0.82 


Ba 


3,6 


HD104979 


n 


/ 


-0.33 


-0.38 


-0.05 


/ 


Ba 


1,9 


HD83548 


/ 


/ 


-0.33 


0.04 


0.37 


0.49 


Ba 


10 


HD65699 


/ 


/ 


-0.3 


-0.45 


-0.15 




Ba 


8 


HR774 


/ 


/ 


-0.3 


0.31 


0.61 


1.03 


Ba 


3,6 


HD116713 


/ 


/ 


-0.29 


0.27 


0.56 


1.62 


Ba 


2 


HD109061 


/ 


/ 


-0.21 


-0.32 


-0.11 


0.6 


Ba 


3 


HD95345 


/ 


/ 


-0.2 


-0.45 


-0.25 




Ba 


8 


HD65966 


/ 


/ 


-0.2 


-0.17 


0.03 


0.49 


Ba 


10 


HD60197 


/ 


/ 


-0.2 


0.21 


0.41 


0.57 


Ba 


3,6 


« Cyg 


/ 


/ 


-0.1 


0.4 


0.5 


0.9 


Ba 


3,11 


HD92626 


/ 


/ 


0.04 


0.68 


0.64 


/ 


Ba 


3 


HD89638 


/ 


/ 


0.1 


0.1 





0.5 


Ba 


3 


^ Cap 


/ 


/ 


0.1 


0.13 


0.03 


1 


Ba 


3,12 


HD121447 


/ 


/ 


0.1 


0.31 


0.21 


0.57 


Ba 


3,6 


HD100012 


/ 


/ 


0.33 


-0.41 


-0.74 


0.28 


Ba 


3 


+42° 2173 


y 


328.3 


/ 


/ 


/ 


/ 


CH 


1 


HD209621 


y 


407.4 


-0.9 


0.11 


1.01 


/ 


CH 


1,13 


+08° 2654A 


y 


571.1 


/ 


/ 


/ 


/ 


CH 


1 


HD5223 


y 


755.2 


/ 


/ 


/ 


/ 


CH 


1 


HD224959 


y 


1273 


-1.6 


0.11 


1.71 


/ 


CH 


1,13 


HD198269 


y 


1295 


-1.4 


-0.49 


0.91 


2.07 


CH 


1,13,14 


HD135148 


y 


1416 


-1.88 


-0.1 


1.78 


/ 


CH 


1,15 
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Table 2 — Continued 



Object 


binaxity 


Period (days) 


[Fe/H] 


[C/H] 


[C/Fe] 


[Ba/Fe] 


class 


Ref. 


HD201626 


y 


1465 


-1.3 


0.01 


1.31 


/ 


CH 


1,13 


HD30443 


y 


2954 


/ 


/ 


/ 


/ 


CH 


1 


HD187861 


/ 


/ 


-1.65 


0.41 


2.06 


/ 


CH 


13 


-38° 2151 


/ 


/ 


-1.4 


-0.39 


1.01 


/ 


CH 


13 


HD189711 


/ 


/ 


-1.15 


0.80 


1.95 


1.80 


CH 


16 


HD42272 


/ 


/ 


-1.10 


1.21 


2.31 


1.65 


CH 


16 


HD 197604 


/ 


/ 


-0.90 


0.85 


1.75 


1.80 


CH 


16 


HD25408 


/ 


/ 


-0.82 


0.70 


1.52 


1.47 


CH 


16 


HD59643 


/ 


/ 


-0.70 


0.58 


1.28 


2.95 


CH 


16 


HD26 


n 


/ 


-0.44 


0.01 


0.45 


/ 


CH 


1,13,17 


HD104340 


/ 


/ 


-1.15 


-0.62 


0.53 


-0.83 


MDBa^ 


3 


HD204613 


y 


878 


-0.35 


0.52 


0.87 


0.56 


SGCH 


17,18,19 


HD89948 


y? 


1153 


-0.27 


0.47 


0.74 


0.83 


SGCH 


17,18,19 


HD122202 


y 


1290 


-0.1 


/ 


/ 


/ 


SGCH 


18 


BD +17° 2537 


y 


1796 


-0.2 


/ 


/ 


/ 


SGCH 


18 


HD202020 


y 


2046 


-0.1 


/ 


/ 


/ 


SGCH 


18 


HD216219 


y 


3871 


-0.32 


0.63 


0.95 
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Note. — 

''FBS: Field Blue Straggler 
''MDBa: Metal-deficient bajrium star 
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(1985); (5) Boffin & Zacs (1994) (6) Smith et al. (1984); (7) Griffin (1991); (8) Barbuy et al. (1992); (9) Luck 



-54- 



(1991); (10) Kovacs (1983); (11) Chromcy ct al. (1969); (12) Smith ct al. (1980); (13) Venture (1992); (14) Lee 
(1974); (15) Carney ct al. (2003); (16) Kipper et al. (1996); (17) Luck & Bond (1982); (18) McClure (1997); (19) 
Smith et al. (1993); (20) Cohen et al. (2004); (21) Preston & Sneden (2001); (22) North et al. (1994); (23) Pereira 
& Junqueira (2003); (24) Sneden et al. (2003) 



[Fe/H] 



Fig. 1. — Sample stars compiled from literature of abundance analysis with high dispersion 
spectroscopy. Filled symbols denote carbon-enhanced, extremely metal-poor (CEMP) stars, 
divided into three evolutionary stages of giants (squares), subgiants (circles), and dwarfs 
(triangles). Open symbols denote CH stars, divided into four groups of CH giants (squares), 
Ba stars (diamonds), subgiant CH stars including both dwarfs and subgiants (circles) and 
field blue stragglers showing strong CH hnes (triangles). Broken line denotes the locus of 
[C/H] = 0. Crosses denote stars without carbon-enhancement, though our compilation is 
not complete for stars of [Fe/H] > —2.5 and for [C/Fe] < 0.5. 
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Fig. 2. — Distributions of carbon abundance, [C/H], among giants (columns enclosed by thick 
lines) and dwarfs + subgiants (shaded columns) of CEMP stars (upper panel) and of CH stars 
(lower panel); in upper panel, CEMP-s + CEMP-nos and CEMP-s giants are enclosed by 
solid and broken lines, respectively. CEMP giants spread in a wide range down to low [C/H] 
values because of dilution by a factor of up to ~ 100, a factor corresponding to the difference 
in the depth of surface convection of giants and dwarfs. For CH giants, the distribution is 
rather narrower since the surface carbon abundance is diluted to be C/0 < 1 because of 
large pristine abundance of oxygen. The spreads in the surface carbon abundance among 
dwarfs and subgiants may be attributed partly to the difference in the carbon enrichment in 
the donor primary stars and partly to the dilution by surface convection. 
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Fig. 3. — Distributions of orbital separations observed from CEMP-s stars (hatched columns) 
and from CH stars (columns enclosed by thick sohd hues) . Dashed and dotted curves depict 
the relative distributions of initial orbital separations, predicted from the binary scenario for 
CEMP-s stars and for CH stars, respectively: a shoulder on the right of each curve comes 
from the contribution of dwarfs. 
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Fig. 4. — Schematic drawings of the hydrogen mixing and progress of hehum-flash driven 
deep mixing (He-FDDM): The hydrogen mixing leading to sphtting of the hehum flash con- 
vective zone to trigger He-FDDM when the helium convection extends through the hydrogen 
containing layer for the first time (top panel): Recurrent hydrogen mixings at small rates 
without splitting the helium flash convective zone, which ends with the strong mixing to 

trigger He-FDDM, and is followed by third dredge-up for stars of mass M > l.SM© (bottom 
n 
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Fig. 5. — Relative abundances between the main s-process elements Ba, or La, and the heavy 
s-process elements, Pb, plotted against the metallicity for CEMP and CH stars. Double 
circles denote those with the orbital periods observed to date; points with names attached 
denote three stars of the shortest periods, which are hkely to have experienced common 
envelope evolution. Broken line denotes the prediction from the radiative ^^C burning model 
by Busso et al. (1999) 
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Fig. 6. — A general picture of evolution of extremely metal-poor stars to the surface carbon- 
enhancement via helium-flash driven deep mixing (He-FDDM) and/or via third dredge-up 
(TDU) on the diagram of the initial mass Mi initial and the pristine metallicity [Fe/H] with 
the CNO abundances in the solar ratio to Fe. The most right-hand side is delineated by 
the upper mass limit, M^^p, to the development of carbon-oxygen, electron-degenerate core, 
taken from Cassisi & Castellani (1993), and hence, giving the lower mass limit to carbon 
ignition under electron non-degenerate condition. 
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Fig. 7. — Enhancement of s-process elements [Ba/Fe], plotted against the enrichment of 
carbon (open symbols) and of nitrogen (filled symbols) for CEMP stars. Squares, circles, 
and triangles denote giants, subgiants, and dwarfs, respectively; diamonds plot three "mixed" 
stars enriched with nitrogen but not with carbon for comparison. 
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Fig. 8. — Evolution to CEMP and CH stars in binary systems on the diagram of the initial 
mass of primary component and the initial separation (the mass of secondary component is 
set at 0.8 Mq). Symbols connected by thin solid lines represent the stellar radii of primary 
stars at the onset of He-FDDM and TDU, taken from the computations of Z = (crosses) 
and [Fe/H] = —3 (boxes) models (Suda et al. 2004, Suda 2006 in preparation) and from the 
computations of more metal-rich populations (diamonds Iben 1975a; Lattanzio 1986; Marigo 
2002). Thick lines denote the orbital separations below which the primary fills its Roche lobe, 
leading to RLO before He-FDDM or TDU starts, and set the minimum initial separations 
for the surface carbon-enhancement. Top three dashed lines denote the separations for which 
the envelope of secondary star is enriched with [C/H] = —3, —2, and —1 through the wind 
accretion; here we assume that the wind is ejected from the primary with the solar carbon 
abundance and that the depth of surface convection of the secondary is 0.35 M0 in mass. 
Four vertical lines denote, from left to right, the lower mass limit to TDU (Lattanzio 1986), 
the upper mass limit to He-FDDM (Fujimoto et al. 2000), and the upper mass limits, M^p, 
to TP-AGB evolution with C + O core for two metallicity of Z = 10^^ and Z = (Cassisi & 
Castellani 1993); and vertical hatched-line indicates the lower mass limit to the hot bottom 
burning in the envelope (Ventura et al. 2001). Models of [Fe/H] = —3 in shaded zones 
become CEMP via RLO or via wind accretion, and models of [Fe/H] = — 1 in hatched zones 
become CH stars via RLO or via wind accretion, respectively. Thin dotted lines denote the 
upper separation limits to become C-rich via RLO. 
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Fig. 9. — Orbital periods vs. surface carbon abundance [C/H] for CEMP stars (top panel) 
and for CH stars (bottom panel). The theoretical predictions from the wind accretion model 
described in §4.2 are plotted under the assumption of Jeans' theorem (solid lines) and of 
constant separations (broken lines) for giants (on the left-hand side) and for dwarfs and sub- 
giants (on the right-hand side) with the primary masses of 1.5 Mq (thin hues) and 3.0 Mq 
(thick lines). The mass of the secondary star is set at 0.8 Mq and the mass in the surface 
convection is taken to be 0.0035 Mq and 0.35 for CEMP dwarf and giants and 0.01 Mq 
and 0.35 Mq for subgiant and giant CH stars, respectively. The carbon abundance in the 
wind is assumed to be [C/H]p = for CEMP stars and [C/H]p = 0.5 for CH stars. Two 
dash-dotted lines in upper panel plot the relations for a wholly convective secondary star of 
mass 0.3 Mq. Symbols denote observational data for CEMP giants (filled squares), subgiants 
(filled circle) and dwarfs (filled triangles) in upper panel and for CH stars (open squares), 
subgiant CH stars (open circles) and field blue stragglers (open diamonds) in lower panel. 
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Fig. 10. — Fractions of CEMP-s and CEMP-nos stars among EMP stars as a function of 
the medium mass Mmd, a parameter of the initial mass function in a log-normal form of 
eq. (36) with the other parameter fixed at Aa^ = 0.33. Thick solid and broken hues plot the 
proportions of CEMP-s and CEMP-nos stars to EMP stars born in binaries, '?/'cEMP-s/V'binary 
and '0CEMP-nos/V'binary, predicted from the binary scenario, respectively: Thin solid lines 
depicts the CEMP-s proportion when we take into account the contribution from low-mass 
stars born as single stars at the equal number with the binaries. Thick dotted line denotes the 
number ratio of CEMP-nos to CEMP-s stars, ^/'cEMP-nos/^CEMP-s- The observed proportion 
of CEMP-s stars (~ 10 ~ 25%) is bounded by two horizontal thin broken lines and two 
shaded areas indicate the parameter ranges of compatible with the observations. Two 
vertical thin dotted lines demarcate the observed bounds of the ratio between CENP-nos 
and CEMP-s stars (~ 1/3 ~ 1). 
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Fig. 11. — The initial mass function of EMP stars on the parameter diagram of the medium 
mass, Mjnd, and the dispersion, Am, in a log- normal form of eq. (36). Thick solid lines depict 
the loci of constant proportions of CEMP-s stars to EMP binaries, '?/'cEMP-s/'?/'binary (= 0.1, 
0.25, and 0.3): thin solid lines denote the loci of constant CEMP-s fractions when the EMP 
stars born as single stars at equal numbers to the binaries are taken into account, and 
hatched area delineates the parameter space of IMFs that give the fraction of CEMP-s stars 
consistent with the observations. Broken lines denote the loci of constant number ratios 
between CEMP-nos and CEMP-s stars, V'CEMP-nos/V'CEMP-s (= 1,0.5,0.33), and shaded 
area defines the parameter space that entails the IMFs compatible with the observations 
both of CEMP-s and CEMP-nos stars. Thin dash-dotted lines denote the loci of constant 
total mass, M[Fe/H]<-2.5, of mother stellar populations (= 10'', 10^ and 10^ Mq) necessary to 
yield such the number of EMP stars, now constituting the Galactic halo, as estimated from 
the recent HES survey under the different assumptions of IMF. Double circle denotes the 
IMF derived by Lucatello et al. (2005b) 
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Fig. 12. — An illustration of relative distributions, m^{m) of stellar masses for the derived 
initial mass function of EMP stars with Mmd = 10 Mq and Am = 0.4; solid and broken 
curves denote the mass distributions of primary components and secondary components, 
respectively. 



